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Abstract 

Stability constants associated with formation of ternary complexes of the type MLA [where M 
= Co(II), Ni(II), Cu(II) and Zn(II), L - 3-Hydroxy-2-methyl-4(H)-pyran-4-one (maltol) or 5- 
Hydroxy-2-(hydroxy methyl)-4(H)pyran~4-one (kojic acid), A = lysine, arginine, histidine, 
aspartic acid and glutamic acid] are determined potcntiometrically at 30°C and 1=0 1M KNO3 
in aqueous medium Ternary complexes of Cu(II) with histidine or aspartic acid suffer loss ot 
stability because of less number of sites available on squaie planar copper compared to 
distorted octahedral copper 

(Keywords: potentiometry / stability constants / maltol / kojic acid / amino acids) 


Introduction 

Specificity leading to selective and preferred interactions is obviously important 
in biological systems. As far as species containing metal ions of 3d series are 
concerned, it is now well established that among mixed ligand complexes those 
composed of the n accepting imidazole group of histidine residue and an O donor are 
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preferably formed 1 ' 3 . Certainly another possibility favouring specificity and selectivity 
are stacking interactions between aromatic residues 4,5 As one may expect that specific 
interactions, for example between nucleotides and amino acids, play a role during the 
first step of biochemical evolution, we felt that it would be of interest to study the 
stability of single ternaiy system consisting of metal ion / ammo acid / oxygen donor 
in order to see if certain ligand combinations are favoured and if any selectivity is 
achieved 

Owing to the biological importance of histidine, we have undertaken the study of 
selected amino acids with transition metals. In continuation of our earlier studies 6,7 , 
we now report the stabilities associated with the metal complexes of selected amino 
acids with bivalent metals in the presence of oxygen donors such as kojic acid or 
maltol. 


Materials and Method 


The ligands kojic acid and maltol were obtained from lancaster (UK). Other 
ligands, lysine, arginine, histidine, aspartic and glutamic acid were obtained from 
Sigma Chemical Company. The metal salts of Co(II), Ni(II), Cu(II) and Zn(II) were of 
AnalaR grade, and their solutions were standardized volumetricaily by titration with 
the sodium salt of EDTA. 8 Potentiometric measurements were carried out at 30°C 
with 0.1 M KN0 3 as background electrolyte in aqueous medium using a control 
dynamics pH meter. The pH meter readings were plotted against a (moles of base 
added per mole of ligand) or m (moles of base added per mole of metal ion). 
Calculations were made with the help of BEST computer programe 9 . 

Results and Discussion 


Binary Complexes : The structures of the ligands and metal complexes are shown 
m Fig. 1. In the binary complexes it is well established that histidine acts as the 
tridentate ligand binding through the carboxylate oxygen, amino nitrogen and 
imidazole nitrogen. But the other tridentate amino acids of the same groups (i.e. arg. 

in ift'l WOU d necessanIy be bidentate with a proton intact (whose pK a is 1 1.98 and 
10.20 for arg- and lys. respectively) on the side chain resulting in the protonated 
comp exes (MHA). Hence, the binary constants refer to values for protonated 
complexes of lys. or arg. The stability constants for normal complexes could not be 
calculated owing to precipitation in the higher pH region. The values obtained are in 
good agreement with the values found in literature 6,7 . 
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Ternary Complexes-. The potentiometric titration curves of ternary complexes, 
MLA where M = Co(II), Ni(II) and Zn(II), L = maltol or kojic acid, A = lys or arg in 
1:1:1 ratio show inflection at m = 1 and m = 2. Hence, the formation of ternary 
complexes are considered in a stepwise manner in which the ligand maltol or kojic 
acid binds first to the metal ion followed by the stepwise addition of amino acids. In 
case of histidine, two inflections obtained one at m = 1 and another at m = 3. The 
ternary complexes of Cu(II) show single inflection atm = 2 (in case of lys. or arg.) or 
m = 3 (in case of histidine). 

The titration curves for glutamic acid or aspartic acid with bivalent metals in 
presence of maltol or kojic acid show inflection at m = 3. The calculations were made 
with BEST computer program 9 and the values thus obtained are given in Table 1. 


The stability of the various ternary complexes and the corresponding binary 
complexes have been compared in terms of the parameter AlogAT. Perusal of AIogAT 
values (Table 1) reveals that the Aiog/f values for ternary complex [maltol-Cu- 
histidine] are more negative (0.72) when compared to the value for [maltol-Cu-lys] or 
[maltol-Cu-arg]. The reason for this observation is due to the fact that in ternary 
complex of maltol or kojic acid and amino acid, the Cu(II) is essentially square planar 
because of the strong ligand field exerted by maltol or kojic acid. In such cases 
histidine would be forced to be bidentate and this would lead to decreased stability 
(i.e. AlogK = -0.72). This sort of problem is not encountered by lys. or arg. as they 
behave bidentate in ternary ?s well as binary complexes. Hence the observed Alog/f 
values (-0.43 or 0.44) have the origin of statistical reasons. Further these values are 
comparable with the values observed for the alanine complex (-0.40). 

Perusal of the AlogAf values also show that aspartic acid is destabilized 
considerably in [asp. + Cu + Maltol] or [asp. + Cu + Kojic acid] ternary complexes 
(i.e. -0.79 or -0.71) when compared to the ternary complexes of alanine (AlogAfc-tOAO 
or -0.45). The reason for destabilization is that asp. in binary complexes binds in a 
tridentate fashion to aquocupric ion, but when it comes to the ternary complexes 
containing strong field ligand like maltol or kojic acid, it becomes bidentate because 
the distorted octahedral Cu(II) ion would be converted to square planar making less 
number of coordinating sites available for asp. 

Comparison of AlogAf values for other metals (Co(II), Ni(II) and Zn(II)) show that 
the ternary complexes of giu. are more stable than the ternary complexes of asp. 
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whether the primary ligand is kojic acid or maltol. The Alog£ values obtained for 
all these metals are consistent with the values obtained on the statistical grounds 
alone. 
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Fig. 1 Ternary complexes of bivalent metals with kojic acid or maltol and selected amino acids 
A - lysine, B - arginine, C = histidine, D = aspartic acid, E = glutamic acid, F = tentative structure 
showing histidine binding as tridentate in binary and bmdentate in ternary complexes 


The Alog£ values obtained for ternary complexes of [asp. + M(II) + KAMA] are 
less negative than values for [asp. + Cu(II) + KAMA]. This shows that all the metal 
ions except Cu(II) have octahedral geometry in which asp. binds to the metal ion 
through three sites as it binds in binary complexes. Hence, in all these metal ions 
except Cu(II), the destablisation due to loss of denticity is not observed. 
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Table 1- logarithms of stability constants of ternary complexes and their corresponding Alog/f values 


[7’=30°C, p 

= 0.1 mol dm' 3 

(KN0 3 ) , aqueous medium] 



Ligands 

Co(II) 

Ni(II) 

Cu(II) 

Zn(II) 


Kojic acid as Primary Ligand (L) 


Alanine 

8.65 (± 0.05) 

10.06 (±0 06) 

14 05 (±0 07) 

9 56 (± 0 02) 


-0 38 

-0 37 

-0 45 

-0 39 

Lysine HC1 

ppt 

9.78 (± 0.06) 

13 73 (±0.08) 

ppt 



-0.40 

-0 44 


Arginine HC1 

ppt 

9 82 ± (0.05) 

13.68 (± 0 03) 

ppt 



-0 39 

-0 41 


Histidine HC1 

12 13 (±0.05) 

12.73 (±0 03) 

15.34 (±0 04) 

1 1 03 (± 0 05) 


-0 42 

-0 41 

-0 68 

-0 43 

Aspartic acid 

9.98 (± 0 05) 

1 1 77 (± 0 06) 

14.19 (±0 03) 

10.17 (± 0.05) 


-0.49 

-0.48 

-0 71 

-0 49 

Glutamic acid 

8 86 (± 0 08) 

10.33 (±0 07) 

13 80 (±0 04) 

9 94 (± 0 06) 


-0 40 

-0.42 

-0 45 

-0.45 


Maltol as Primary Ligand (L) 


Alanine 

9 28 (± 0 08) 

10 77 (±0 07) 

15 14 (±0 05) 

10.26 (±0 06) 


-0 35 

-0.30 

-0.40 

-0 32 

Lysine HC1 

ppt 

10.45 (±0 08) 

14.78 (±0 07) 

ppt 



-0 37 

-0.43 

- 

Argmme HC1 

ppt 

10 46 (±0 07) 

14 69 (±0 06) 

ppt 



-0 39 

-0 44 


Histidine HC1 

12 72 (± 0.07) 

13 41 (± 0.08) 

16 34 (±0 08) 

1 1 67 (± 0 07) 


-0.43 

-0 37 

-0 72 

-0 42 

Aspartic acid 

10 60 (± 0 06) 

12 48 (±0 08) 

15 15 (±0 03) 

10 84 (±0 06) 


-0.47 

-0 41 

-0 79 

-0 45 

Glutamic acid 

9 46 (±0 08) 

11 00 (±0 06) 

14.82 (±0 04) 

10.61 (±0 05) 


-0.40 

-0 39 

-0 47 

-0 41 


Sigma fit values are given in the parentheses. 
AlogK = log Ca - [log C. + log C'] 
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Abstract 

Mixed ligand complexes of the type [ML(Aa)], where M = Co(II), Ni(II), Cu(II) and Zn(II), L 
= kojic acid or maltol, Aa = amino acids, viz glycinate, alaninate, a-aminobutyrate, norvalme, 
norieucine, vahnate, leucinate or isoleucmate have been studied by potentiometnc pH titrations. 
The potentiometric measurements reveal a slightly higher formation tendency for the system 
with leucmale as compared to those with alaninate This increase in stabilization is attributed to 
an intramolecular hydrophobic ligand-ligand interaction between n electron cloud of kojic acid 
or maltol and the isopropyl group of leucine. 

(Keywords : potent! ometry / ternary complexes / formation constants / maltol / kojic acid / 
amino acids) 


Introduction 

Hydrophobic interactions are known to occur among other noncovalent 
interactions in biomolecules, 1 they contribute to the formation of distinct structural 
features 2 Such hydrophobic interactions occur between two aliphatic groups or 
between aliphatic and aromatic moieties; the ring stacking between two aromatic 
systems is also often listed within this category 3,4 

In continuation of our earlier studies 5,0 on the hydrophobic studies, we now report 
the stability constants associated with ternary complexes containing kojic acid or 
maltol as primary ligand and bivalent metals in presence of amino acids with aliphatic 
side chains such as glycine, alanine, cx-ami nobutyric acid, norvaline, norieucine, 
valine, leucine or isoleucinu. These amino acids were selected such that the side chain 
varied in length and volume from a methyl group as in alanine to n-butyl group as in 
norieucine or an isopropyl group as in leucine. 
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It is oovious that the methyl side chain in alanine is too short to allow any 
interaction with the aromatic rings within the square planar Cu(II) complexes and for 
the octahedral metal complexes, it is expected to be insignificant if it does occur at all 
The reason why alanine and not glycine was selected as the basis for the comparison 
is that Martin 7 has shown that the stability of the glycine complexes is often some 
what exceptional. The alanine complexes are therefore a better basis for comparison. 

Materials and Method 

The ligands, maltol, kojic acid were obtained from Lancaster (U.K). The amino 
acids were obtained from Sigma Chemical Company. The metal salts of Co(II), Ni(II), 
Cu(II), and Zn(II) were of AnalaR grade and their solutions were standardized 
volumetrically by titration with the sodium salt of EDTA. 8 Potentiometric 
measurements were carried out at 30 C in aqueous medium (with 0.1 M KN0 3 as 
background electrolyte) using control dynamics pH meter. The pH meter readings 
were plotted against a (moles of base added per mole of ligand) or m (moles of base 
added per mole of metal ion). Calculations were carried out with the help of computer 
program - BEST. 9 


Results and Discussion 

Binary stability constants pertaining to the interaction of amino acids, kojic acid 
and maltol with bivalent metals were redetermined to have a better comparison 
between binary and ternary complexes under the same conditions and the values thus 
obtained are in good agreement with the literature values 5,6 . 

The potentiometric titration curves of [ML(Aa)] in 1:1:1 ratio where L = maltol or 
kojic acid; Aa = amino acids under study; M = Co(II), Ni(II), and Zn(ll), show a 
inflection at m =1 followed by another at m = 2. The ternary curve coincides with 
binary curve upto m = 1 where 1:1 [ML] complexes are formed. This indicates that 
only 1 : 1 [M-L] complexes are formed in the region m = 0 to m =1 and amino acid does 
not coordinate with M(II) in this region. The divergence of mixed ligand curve for the 
binary [M-L] curve beyond m — 1 suggests the coordination of amino acid as a 
secondary' ligand. Hence, the formation of ternary complexes can be considered to 
occur in stepwise equilibria. The ternary constants thus obtained are presented in 
Table 1. In case of Cu(II) complexes in 1:1:1 ratio, the titration curve shows an 
inflection at m — 2 indicating the simultaneous complexation. 
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Table 1- logarithms of stability constants of ternary complexes 
[T = 30°C, ji = 0 1 mol dm' 3 (KNQ 3 \ aqueous medium] 


Ligands 

Co(II) 

Ni(II) 

Cu(II) 

Zn(Il) 



Kojic acid as primary ligand (L) 


Glycine 

8 94 

10.37 

14 12 

9 64 

a-aminobutyric acid 

8 60 

9.79 

14 20 

9 32 

Norvalme 

8.63 

10.39 

14.42 

9 53 

Norleucine 

8 70 

10 17 

14 58 

9 42 

Valme 

8 68 

10 36 

14 36 

9 50 

Leucine 

8 70 

10 23 

14 33 

9 60 

Isoleucme 

8 65 

10 19 

14.29 

9 50 



Maltol as primary ligand (L) 


Glycme 

9 57 

11 13 

15 18 

1032 

a-aminobutyric acid 

9 24 

10.51 

1531 

10.04 

Norvahne 

9 27 

11 11 

1551 

10 27 

Norleucine 

9 35 

10 89 

15 68 

10 15 

Valine 

9 35 

11 10 

15.49 

10.28 

Leucine 

9.35 

10.97 

15 46 

10 33 

Isoleucine 

9.31 

10.93 

15 40 

10.24 


Sigma fit value are in the range ±0 05 to 0 08 log units 

Table 2 lists the Alog K values along with AAlog K values. The parameters AlogAT 
and AAlogAT are used to characterise the stability of ternary complexes formed. 
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Negative Alog£ values suggest that the binary complexes are more stable than the 
ternary complexes, whereas a positive MogK suggests greater stability of the ternary 
complexes. The magnitude of Alog K values depends on the denticities of both the 
ligands and coordination number of metal ion. In the present study kojic acid or maltol 
and all amino acids are listed as bidentate ligands. If both ligands in ternary 
complexes are bidentate and geometry is octahedral, then the statistical factor would 
be 5/12 and its corresponding log value shall be -0.4. Thus if the obtained AlogK 
values are positive or less negative than the statistical value (i.e. -0.4) then ternary 
complexes are more stable. However, it should be noted heie that negative values of 
AlogK do not preclude the formation of ternary complexes in solution. 

Table 2- Comparison of Alog£ and AAIog/f values for ternary complexes 


[r=30°c,u= 

0 1 mol dm' 3 

(KN0 3 ), aqueous medium] 



Secondary Ligands (A) 

Co(II) 

Ni(ll) 

Cu(II) 

Zn(II) 


Kojic acid as primary ligand (L) 


Glycme 

-0 43 

-0 48 

-0 46 

-0 43 

a-ammobutyric acid 

-0.36 

-0 33 

-0.39 

-0 35 


(+0 02) 

(+0 04) 

(+0.06) 

(+0.04) 

Norvahne 

-0.32 

-0 29 

-0.33 

-0.33 


(+0.06) 

(+0 08) 

(+0.12) 

(+0 06) 

Norleucine 

-0 28 

-0 24 

-0.30 

-0 30 


(+0 10) 

(+0.13) 

(+0 15) 

(+0.09) 

Valine 

-0 32 

-0.28 

-0.33 

-0.29 


(+0.06) 

(+0 09) 

(+0.12) 

(+0 10) 

Leucine 

-0 26 

-0.22 

-0 28 

-0 26 


(+0.12) 

(+0 15) 

(+0.17) 

(+0.13) 

Isoleucine 

-0 35 

-0.32 

-0.37 

-0.34 


(+0.03) 

(+0.05) 

(+0 08) 

(+0.05) 

Table 2 Contd... — 
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Table 2 Contd. . 

Maltol as primary ligand (L) 


Glycine 

-0 40 

-0 36 

-0 44 

-0 38 

a-aminobutyric acid 

-0 32 

-0 25 

-0 32 

-0.26 


(+0 03) 

(+0.05) 

(+0.08) 

(+0.06) 

Norvahne 

*0.28 

*0 21 

-0.28 

-0.22 


(+0.07) 

(+0.09) 

(+0 12) 

(+0 10) 

Noi leucine 

-0 23 

-0 16 

-0 24 

-0.20 


(+0 12) 

(+0 14) 

(+0 16) 

(+0 12) 

Valine 

-0 25 

-0 18 

-0 24 

-0 14 


(+0.10) 

(+0 12) 

(+0 16) 

(+0 18) 

Leucine 

*0 21 

-0 12 

-0 19 

-0 16 


(+0 14) 

(+0 18) 

(+0 21) 

(+0 16) 

Isoleucine 

*0 29 

-0.22 

-0 30 

-0 23 


(+0 06) 

(+0 08) 

(+0.10) 

(+0 09) 


Alog K = log - [log K^ a + log ] The AAiogA values are presented in paranthesis. 

AAlogA = [AIogA, - Alog/y, where Alog A, is for the reaction ( M + L + amino acids except alanine) 
and AiogA' 2 is for the reaction (M + L + alanine). 

Alanine is taken as a reference for zero based hydrophobic interactions 

The determined stability constants of the ternary complexes i.e., especially the 
values of AAlogAT, give hints into the direction of an intramolecular hydrophobic 
ligand-ligand interaction in M(L)(Aa) complexes. A tentative structure of such a 
ternary complex is shown in Fig. 1. In this structure the isopropyl group of the amino 
acid is located above the n electron cloud of kojic acid / maltol, and the two ligands 
are linked together by the metal ion; in other words, the metal ion promotes the 
hydrophobic interactions. 
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Extent of the intramolecular ligand-ligand interactions 

Based on the results, it can readily be understood why the M(L)(Aa) complexes 
show enhanced stability: the isopropyl group of amino acid is ideally suited for 
hydrophobic ligand-ligand interactions The occurence of complexes with a structure 
similar to other shown in Fig. 1, which is responsible for the slight increase in stability 
does not imply that all the M(L)(Aa) ternary species exist in this closed form. 
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€00“ 
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H,N— C -H 
$ I 
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ch 3 
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C0G~ 
+ I 

H-N-C-H 

s 

CCHi 3 

CH 3 

E 
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C00~ 
+ i 

h 3 n-c-h 


H3C 


CM 

X ''CH3 


Hydrophobic interactions 


C00~ 

i 

h 3 n~c-h 

i 

ch 2 

1 

CH 

h 3 c^ Nh ch 3 
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COO* 
+ I 

HnN-C-H 
4 I 


H3C 


CH3 



7 “ “ u-ammoDutync acid, D = norvaline, E = norleucme, F = valme G = 

leucine, H - isoieucine. 1 = tentative structure showing the intramolecular hydrophobic 
interactions in the ternary complex of [maltol-M-leucme] 

i„J"Jf' U ! i 0 ",’l’ | e h re is “ rtain| y “ intramolecular, and therefore eoneentration 
ndependent, equilibrium between an 'open' and a closed form as shown below. 
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The dimensionless constant K\ of this intramolecular equilibirum is defined by the 
equation 1. 


K x = 10 AAlo §^-i 


( 1 ) 


The other parameter AAlogA - is given by relationship 

AAlog£ = (Alogtfi - Alog£ 2 ) (2) 

where AlogA] is for the reaction [M + L + Aa] and AiogA2 relates to the reaction 
[M + L + alanine]. 

The percentage of stacked isomers in the systems are also calculated with the 
relationship 3: 

%(MLA) hj = [100 * / ('I + if,)] (3) 

The values of K\ and % (MLA) H i for these systems are presented in Table 3. The 
perusal of parameters associated with hydrophobic interactions reveal that 
stabilization in ternary complexes with respect to amino acids follow the increasing 
order : 

Straight chain: 

a-aminobutyric acid <Norvaline <Norleucine 
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Branched chain 

Isoleucine < Valine < Leucine 

It is obvious from the above resuits that the extent of the intra-molecular 
hydrophobic interactions in these ternary complexes are dependent on a) length and 
b) possibly also on the bulkiness of the aliphatic side chain of the amino acid. 

With increase in chain length (viz., as propyl and butyl groups in norvaline or 
norleucine) there is better overlap between two ligands which results in maximum of 
hydrophobic interaction as reflected by their parameters. 

Comparison of % stacked isomers formed (Table 3) for leucine and norvaline 
system reveal that leucine (having branched side chain) forms better stacked isomers 
than the norvaline, though both ligands have the same number of carbon atoms. This 
is an evidence that the bulkiness of aliphatic side chain of the amino acid play a role in 
enhancing the hydrophobic interactions in these ternary complexes. 

Table 3- Parameters showing the extent of intramolecular hydrophobic interaction in ternary complexes 


[T = 30°C, p = 0. 1 mol dm' 3 (KNO } ), aqueous medium] 


Secondary 
Ligands (A) 

Parameters 

Co(il) 

Ni(JI) 

Cu(II) 

Zn(II) 


Kojic acid as pi imary ligand (L) 



a-ammobutyric 

Kt 

0 047 

0 096 

0 148 

0 096 

acid 







%(MLA) HI 

4 48 

8 75 

12 89 

8.75 

Norvaline 

K 1 

0 148 

0 202 

0318 

0 148 


%(MLA) h1 

12 89 

16 80 

24 13 

12.89 

Norleucine 

K 

0 258 

0 348 

0412 

0 230 


%(MLA) HI 

20 50 

25 82 

29 17 

18 70 

Valine 

*1 

0 148 

0 230 

0318 

0 258 


%(MLA) Hi 

12.89 

18 70 

24.13 

20.50 

Leucine 

Ki 

0318 

0412 

0 479 

0 348 


%(MLA)st 

24.13 

29 178 

32.38 

25.82 


Table 3 Contd 
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Table 3 Contd 


Isoleucme 

Ki 

0 071 

0 122 

0 202 

0 122 


%(MLA) hi 

6 67 

10 87 

16 80 

10 87 



Maltol as 

primary ligand (L) 



a-aminobutync 

acid 


0 071 

0 122 

0 202 

0 148 


%(MLA) hi 

6 67 

10 87 

16 80 

12 89 

Norvahne 

K , 

0 174 

0 230 

03182 

0 258 


%(MLA) hi 

14 82 

18 70 

24 13 

20 50 

Norleucine 

Ki 

0318 

0 30 

0 445 

0318 


%(MLA) hi 

24 13 

27 54 

30 79 

24 13 

Valine 

Ki 

0 258 

0318 

0.445 

0513 


%(MLA) H! 

20 50 

24 13 

30 79 

33 91 

Leucine 

K t 

0 380 

0513 

0.622 

0.445 


%(mla) hi 

27 54 

33 91 

38 35 

30.79 

Isoleucine 

Ki 

0 148 

0 202 

0 258 

0 230 


%(MLA) hi 

12.89 

16 80 

20 50 

18 70 


Further, the values of % stacked isomers of isoluecine complexes are more or less 
same as that of values of a-aminobutyric acid complexes, indicating that though 
isoleucine has branched side chain it behaves as a-aminobutyric acid in ternary 
complexes (i.e. it appears to react as an ethyl substituted a-aminobutyric acid). Hence 
it is concluded that in hydropnobic interactions, steric restriction does exist. The order 
with respect to primary ligands (maltol > kojic acid) is in line with basicity of the 
ligands. Again order with respect to metals is in confirmity with Irving Williams 
natural order of stabilities 10 . 
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Abstract 

A few Mn(II), Cr(III) and Fe(III) complexes have been prepared from 4,4'-diaminostilbene- 
2,2'-disuIphonate (DAS) with • anillin (L = DASV) and o-vamllin (L' = DASOV) Schiff bases 
These complexes confine to stoichiometry of the type [Mn 2 L 2 Cl 4 (H 2 0) 4 ], [Fe 2 L 2 C1 6 (H 2 0) 2 ], 

[Cr 2 L 2 CI 6 (H 2 0) 2 ] and [Mn L' 2 (H 2 0) 4 ], [Fe 2 L' 2 Cl 2 (H 2 0) 2 ], [Cr 2 L' 2 (H 2 0) 2 ] The complexes 
have been characterized by elemental, conductance, magnetic and spectral data. The infrared 
spectra of the complexes suggest that nitrogen and oxygen are involved in the bonding m case 
of DASOV ligand. However in DASV ligand the OH groups are not involved in bond 
formation The magnetic and electronic spectra indicate that Mn(II), Cr(IIi) and Fe(III) 
complexes have octahedral configuration 

(Keywords IR/UV/stilbenediamine Schiff bases/Mn(II), Cr(III)/Fe(III)) 

Introduction 

There have been a number of studies on metal complexes of Schiff bases formed 
from salicylaldehyde, substituted salicylaldehydes with various aromatic amines 1 ” 5 . 
4 s 4'-diaminostilbene-2,2'-disulphonate and its derivatives with vanillin and o-vanillin 
constitute a very important class of compounds because they are biologically 
important. Substituted stilbene diamine derivatives have been found to posses antiviral 
activity 6 . The literature 7 study shows that there appears no reports on corresponding 
chelates of Mn(II), Cr(III) and Fe(III) with ligands of stilbene diamines. Consequently 
the present work deals with the study of nature of metal-ligand bonding by 
spectrochemical and magnetochemical data with the following ligands (DASV and 
DASOV). 
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(DASV) (R = S0 3 Na) 



(DASOV) (R = S0 3 Na) 

Materials and Method 

The chemicals used in the preparation of ligands were of reagent grade. The 
following metal salts (MnCI,.4H 2 0, FeCi 3 .6H 2 0, CrCl 3 .6F5,0) were used foi 

preparation of metal complexes. 

Preparation of ligands : In a typical preparation to methanol solution of (0 04 
mole) vanillin was added 0.02 mole of 4,4'-diaminosti!bene-2,2'-disu!phonate and tne 
reaction mixture was refluxed for about 2 h. The Schiff base (DASV) separated aftei 
cooling was filtered and recrystallised from ethanol The same procedme was 
followed for o-vanillin Schiff base (DASOV) 

Preparation of metal complexes : To an ethanolic solution ot 0 01 mole of DASV 
/ DASOV was added 0.01 mole of MetalGI / 111) chloride and the reaction mixture 
was refluxed for about 1 h and then sodium acetate (2g) was added to the reaction 
mixture and the refluxing was continued for about 3 h. The complex so obtained was 
filtered, washed with water and then with alcohol and finally dried in vacuum over 
fused calcium chloride. 
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Analysis : The complexes were analyzed for their metal, chloride and nitrogen 
contents by standard methods 8 The results of the elemental analysis are shown in 
Table 1. 

Physical Measurements * The molar conductance measurements of the complexes 
m DMF at a concentration ~10’ 3 M were taken with an ELICO conductivity bridge 
having a cell of cell constant 0.829 cm" 1 manufactured by Electronic Industrial 
Instruments Co., Hyderabad. The magnetic measurements were taken at room 
temperature on a Gouy balance with glass tube calibrated for diamagnetism. The 
electronic spectra were obtained with a Beckman DB-2 recording spectrometer using 
1 cm glass cell in the region 350-1200 nm. The weighed samples were dissolved in 
DMF and diluted to the required concentration ~10‘ 4 M. The IR spectra of the ligands 
and the complexes in KBr pellets were recorded in a Perkin Elmer 621 
spectrophotometer in the range 4000 - 400 cm’ 1 . The Sow region IR spectra in the 
range of 450 - 250 cm’ 1 were recorded in polyethelene pellets. 

Results and Discussion 

The Mn(II) and Fe(III) complexes were reddish brown, where as those of Cr(III) 
were reddish in colour. All these complexes have 1:1 stoichiometry (Table 1). They 
were insoluble in common organic solvents. The molecular weight of these complexes 
could not be determined because of their insolubility in nitrobenzene. However, they 
were soluble to a limited extent in DMF, DMSO and pyridine. The molar conductance 
values in DMF were too low to account for any dissociation. Hence, they can be 
regarded as non-electrolytes 9 . * 

Magnetic data : The magnetic moments obtained at room temperature for these 
complexes are listed in Table 1. 

The magnetic moments for manganese(Il) complexes are in the range of 5.35 to 
5.85 B M. suggesting the formation of spin-free complexes. This closely agrees with 
the expected value of about 5 90 B.M for Mn(II) complexes 10 . 

The chromium(III) complexes exhibit magnetic moments in the range of 3.55- 
3.87 B.M agreeable to spin only value of chromium(III) complexes with octahedral 
configuration, viz., 3.88 B.M 11 . 

The iron(III) complexes are expected 12 to possess magnetic moments very close to 
spin only value of 5.92 B.M. The magnetic moments obtained for these iron(III) 
complexes were in the range of 5.89-5.90 B.M. All these complexes have lower than 
the spin-only values because of polymerization. 



Table I— Elemental analysis, colour, conductance, magnetic and electronic spectral data for DASV and DASOV and their met. 
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Electronic spectra : The observed band maxima for the complexes are listed in 
Table 1. The electronic spectra of manganese(II) complexes show high intensity band 
maxima observed in the 24000-22480 cm' 1 region which has been regarded as ligand 
metal charge transfer band. The remaining d-d bands are observed in the region 
20000-19000 and 13000 cm' 1 l3 . The Fe(III) complexes show a high intensity band 
around 25000 cm' 1 due to the ligand metal charge transfer, this obscures the other d-d 
transitions in the visible region 14 . It has been observed in the diffuse reflectance 
spectra that iron(III) complexes exhibit two bands in the 12500-10640 and 20000- 
16670 cm' 1 regions attributable to 6 X, g 4 7’ |g (G) and 6 A lg -> 4 A lg (G ) transitions 

respectively 15 . In these complexes we observe weak and broad band maxima in the 
regions 12500-11750 cm' 1 and 16140-16000 cm" 1 and these can be attributed to the 
6 A ls — > 4 r, g (G) and 6 A lg ~> 4 A lg (G) transitions. 

Octahedral chromium(III) generally exhibits two spin allowed bands in the visible 
region 16 . The low energy transition 4 A 2g -> 4 J 2 g occurs in the 17000-19000 cm' 1 
regions and another band due to 4 A 2g —> 4 T\ g (F) transition is present in the region 
27000-30000 cm' 1 . In these complexes of chromium(III), we observe two bands in the 
following regions 25000-21720 cm' 1 and 18400-19230 cm' 1 . These have been 
attributed respectively to 4 A 2g -> 4 7’i k (F) and 4 A 2g -> 4 T 2g transitions. The high intensity 
of the first band may be due to overlaying of the charge transfer band. 

Infrared spectra : Important infrared frequencies of ligands and their complexes 
are set out in Table 2. It is shown 17 that free hydroxy group gives rise to stretching 
band in the region 3600-3500 cm' 1 . This band shifts to lower frequency on hydrogen 
bond formation. The intermolecular hydrogen bonding shifts it further to the lower 
frequency and causes broadening. The intramolecular hydrogen bonding further shifts 
it to the lower frequency. Hence, the band becomes weak and broad taking these 
observation into consideration and an analogy with salicylideneaniiine, the weak 
broad band observed around 2650 cm' 1 is assigned to H-bonded -OH. 

Considering the previous assignments 18 the strong band observed around 1620 
cm' 1 for Schiff base has been assign to v(C=N) vibration. In all these complexes this 
band appears at lower frequency region 1620-1600 cm' 1 indicating that the coordinate 
bond is formed between the nitrogen of the azomethine group and metal ion. 

There are three bands of high intensity in the range 1600-1550 cm' 1 which are 
attributed to the aromatic ring vibrations 19 . These bands show remarkable changes in 
frequencies on going from the ligand to metal complexes. 



Table 2- Important infrared frequencies (cm’ 1 ) of DAS V and DASOV and their metal complexes 
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* Intramolecular -H-bonded -OH 
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A strong intensity band observed in the region 1280 cm" 1 of Schiff base ligand 
(DASV) has been assigned to the phenolic v(C-O) vibration. This band does not show 
any perturbation m the complexes indicating non involvement of OH groups in the 
bond formation. These observations emphasize that the ligand (DASV) has not 
undergone deprotonation while reacting with metal(II/III) chlorides. Analytical data 
also support this view. The asymmetric v(C-Q) vibrations of the -OCH 3 group of the 
ligands at 1320 cm' 1 remains inert to the complex formation indicating non 
involvement of OCH 3 group of ligand in coordination to the metal ion. 


However a strong intensity band observed in the region 1280 cm" 1 of Schiff base 
ligand (DASOV) has been assigned to phenolic v(C-O) vibration. In view of the 
previous 20 assignments, the H-bonded -OH band around 2700 cm* 1 has been assigned 
to the intramolecular -K bonded -OH m this ligand. Further the assignment of the 
band is confirmed by comparing the spectra of the benziledene aniline and anisiiidene 
aniline with that of DASOV ligand. The band due to intramolecular -H bonded -OH is 
absent in the model compounds. This band disappears in the spectra of complexes 
indicating reaction through oxygen via deprotonation. Further, the band at 1280 cm’ 1 
attributed to v(C-O) has been shifted to higher frequency and is found in the region 
1340-1305 cm’ 1 for these complexes. This confirms the involvement of phenolic -OH 
in the complex formation. The complexes display a broad intense band centered 
around 3400 cm' 1 attributed to water molecules. No loss in weight was observed even 
after drying the complexes at 105°C for 2 h. The infrared spectra of the dried 
complexes also exhibit bioad intense band around 3400 cm' 3 and -830 cm' 1 . 


The high intensity bands appearing in the region often interfere with M-L bands 
Thus assignments of the bands to various v(M-N), v(M-O), v(M-Ci) vibrations m this 
region become much complicated. However, the assignments of the bands to vailous 
modes have been made by comparing the spectra of the complexes with those of 
ligands Taking into consideration of previous reports 25 the non-ligands bands in the 
region 543-520 cm’ 1 , 460-440 cm’ 5 and 348-318 cm’ 3 have been assigned to v(M-N), 
v(M-O) and v(M-Cl) vibrations respectively. 


The elemental analysis, magnetic study, electronic and IR spectral observations 
project the following structures of coordination number six for these complexes. 
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(R = S0 3 Na,X = H 2 0) 
[Mn 2 (DAS0V)2(H 2 0)4] 



[Mn 2 (DA SV) 2 Cl4(H 2 0)4] 



M = Fe(III), Cr(ni), X = H 2 0 

[Cr 2 (DAS0V),Cl 2 (H 2 0) 2 ] [Fe 2 (DAS0V) 2 CI 2 (H 2 0) 2 ] 



SYNTHESIS & CHARACTERIZATION OF Mn(II), Cr(III) & Fe(III) COMPLEXES 243 



( R = S0 3 Na, X = H 2 0 ) 
[Fe 2 (DASV) 2 Cl 6 (H 2 0) 2 ] 
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Abstract 

A computer assisted pH metric investigation has been made on the speciation of Ca(II) and 
Mg(II) complexes of glutamine and succinic acid in varying concentrations (0-60% v/v) of 
acetomtrile-water solutions at an ionic strength of 0 16 mol dm" 3 (NaCl) and 303 K 
temperature The predominant species detected for Ca(II) and Mg(II) complexes of glutamine 

2 

are ML 2 , ML 2 H + and of succinic acid are ML^' , MLH + and ML 2 H The variation in values of 

stability constants with the change in dielectric constant of medium is attributed to the 
electrostatic and non-electrostatic interactions 

(Keywords speciation study/complex equilibria/acetonitrile/L-glutamine/ succinic acid/ 
speciation) 


Introduction 

The physiological activities of Gin, Sue, Ca(II) and Mg(II) are associated with the 
metabolic processes in liver; hence, their speciation has been studied in varying 
compositions of acetonitrile-water mixtures. The protonation equilibria of Gin and 
Sue in acetonitrile-water have been reported earlier 1 . 

Materials and Method 

Solutions of L-glutamine, succinic acid, calcium(II) and of magnesium(II) 
chlorides (E. Merck, Germany) were prepared in triple distilled water. 99.5% pure 
acetonitrile (E.Merck) was used. The strength of alkali and the acid present in the 
metal ion solutions was determined using the Gran plot method 2 . The titrations were 
carried out in the medium containing varying concentrations (10-60% v/v) of 
acetonitrile in water maintaining an ionic strength of 0.16 mol dm' 3 with NaCl at 
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303.0 ± 0.1 K using a Systronic (model 335) pH meter of 0.01 readability. In each of 
the titrations, the titrand consisted of mineral acid of approximately 1 mmol in a total 
volume 50 cm 3 . Titrations were carried out with 0.4 mol dm' 3 sodium hydroxide. 
Other experimental details are given elsewhere 3 . The best fit chemical models for each 
system investigated were arrived at using a non-linear least squares computer program 
MINIQUAD75 4 


Results and Discussion 

Refinement of species : The final models refined for Gin and Sue complexes 
contained ML 2 and ML 2 H + and ML 2 , MLHf and ML 2 H\ respectively (Tables 1 and 
2). A very low standard deviation in log (3 values indicates the precision of the 
parameters. The deviation of residuals from Gaussian distribution is reflected by the 
statistical parameters, skewness and kurtosis Skewness indicates whether the 
residuals are concentrated to the left or to the right of the mean. Kurtosis indicates the 
peaked ness of the residual curve. For an ideal normal distribution, the values of 
kurtosis and skewness should be three and zero, respectively. The kurtosis values of 
most of the systems indicate that the residuals form platykurtic pattern. The values of 
skewness evince that the residuals of many systems form a part of normal distribution. 
The sufficiency of the model is further evident from the low crystallographic R values. 

Effect of errors in influential parameters like the concentrations of ingredients on 
the stability constants is studied and the data show that the order of affecting the 
magnitudes of stability constant is due to errors in the concentrations of alkali> 
ligand> acid> metal. 

Solute-solvent interactions: The linear variation of logK values of Ca(II) and 
Mg(II) complexes of L-glutamine and succinic acid with 1 ID indicates that 
electrostatic forces are dominating the equilibrium process under the present 1 
experimental conditions. The decrease in the stability constants of Mg(II) complexes 
of succinic acid at 60% may be due to the competition of solvent molecules with the 
ligand for coordination position of metal ion in higher percent of acetonitrile-water 
mixture. Since complex formation can be viewed as a competition between the pure 
and solvated forms of ligand and the metal ion, the solute-solvent interactions, relative 
thermodynamic stabilities and kinetic labilities are also expected to play an important 
role. Different types of electrostatic forces dominate in different ranges of the 
composition of acetonitrile-water mixtures. The deviation from linearity in Ca(II)-Suc 
indicates the blend of both electrostatic and non-electrostatic interactions. 
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Table 1- Best fit chemical models of Ca(II) and Mg(II) complexes of L-glutamine in acetonitrile-water 
mixtures. 

No. of titrations in each percentage = 4 


% v/v 

AN 

lOgPmlh(SD) 

120 121 

NP 

Ucor, X 

10 7 

Skewness 

Kurtosis 

x 2 

R 



CafllV Glutamine (dH ranee '8 5-11 0) 



00 0 

3 54(5) 

13.46(7) 

28 

2 577 

-0 17 

78 

35 41 

0.015 

10 0 

3 46(7) 

13 91(1) 

27 

2 874 

-1 77 

10 04 

31 30 

0012 

20.0 

3 80(1) 

13 86(9) 

28 

3 790 

-2 94 

17 96 

76.00 

0013 

30.0 

3 72(6) 

13 90(8) 

28 

2 824 

-3 14 

18 30 

60 95 

0.011 

40 0 

3 89(7) 

13.85(2) 

27 

5 192 

-2 14 

9 04 

28 33 

0015 

50 0 

401(1) 

14.26(1) 

28 

2 902 

3 52 

-21.59 

96 38 

0012 

60 0 

3 96(9) 

14.41(8) 

28 

2 167 

-3 01 

18 40 

46 29 

0010 



Medn-Glutamine (pH range 8 3 - 10 6) 



00 0 

3 .40( 5) 

12.92(9) 

28 

8 185 

-0 11 

64 

67 38 

0 025 

10 0 

3.42(2) 

13 55(1) 

28 

1 826 

-1 48 

14 88 

78 86 

0 009 

20 0 

3 61(4) 

13.62(3) 

27 

0 503 

-0 61 

6.46 

30 70 

0 004 

30 0 

3.23(8) 

13.72(4) 

26 

0.369 

-0 24 

5.36 

15 13 

0 004 

40 0 

3 73(7) 

13 78(8) 

28 

1 348 

-1 90 

9 26 

42 10 

0 008 

50 0 

4.07(5) 

14.10(6) 

27 

1 154 

-1 92 

8.53 

40 78 

0 007 

60 0 

3.36(2) 

14.21(3) 

27 

0 296 

1.94 

14 07 

73 17 

0 003 
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Table 2- Best fit chemical models of Ca(II) and Mg(II) complexes of succinic acid in acetonitrile -water 
mixtures 

No of titrations in each percentage = 4 


% 

v/v 

AN 

120 

10gPm!h(SD) 

111 

121 

NP 

UcorrX 

10 7 

Skewness 

Kurtosis 

X 2 

R 




CadlV Succinic acid 

(dH ranee 

3 0-9 0) 




00 0 

3 40 

6 36(5) 

9 61(8) 

621 61(8) 

1 626 

-0 06 

6 1 

62 78 

0 034 

100 

Rejected 

Rejected 

9 74(9) 

61 

3 200 

-2 22 

13 22 

57 66 

0 008 

20 0 

4 00(7) 

7 18(2) 

9 99(7) 

62 

1 406 

-0 57 

4 37 

45 05 

0 006 

30 0 

4 43(4) 

7 64(4) 

10 51(5) 

57 

0 948 

-0 97 

6 87 

43 78 

0 005 

40 0 

3 65(7) 

7 71(3) 

10 27(7) 

58 

1 531 

-0 31 

5 84 

57 45 

0 006 

50 0 

3 61(6) 

Rejected 

10 66(6) 

59 

1 025 

-1 17 

6 65 

68 83 

0 005 

60 0 

4 24(4) 

8 17(3) 

1 1 00(4) 

60 

0 586 

-1 23 

6 74 

43 82 

0 004 




MsflIVSuccmic acid (dH ranee 

3 0 - 9 0) 




00 0 

3 66(4) 

7 16(1) 

961 

62 

4 498 

-3 76 

23 48 

109 7 

0016 

100 

Rejected 

Rejected 

9.64(6) 

61 

8 136 

-1 46 

8 95 

1186 

0 014 

20 0 

4 23(5) 

7 49(1) 

10 06(6) 

62 

1 337 

-0 82 

5 14 

43 38 

0 005 

30 0 

3 96(5) 

Rejected 

10 11(1) 

57 

2316 

-1 50 

8 06 

59.22 

0 008 

40 0 

4 22(3) 

7 62(3) 

10 41(4) 

59 

0 576 

-0.75 

5 98 

57 80 

0 004 

50.0 

4.21 (2) 

8.06(6) 

10.78(2) 

60 

0 268 

-1 75 

10.20 

56 80 

0 002 

60 0 

3 67(3) 

Rejected 

10 76(4) 

61 

0 562 

-1 77 

32 23 

85.20 

0.003 
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Distribution diagrams : Different forms of Gin are LH 2 , LH, L" in the pH regions 

1.54.5, 1. 5-9.5 and 8.0-11.0, respectively 1 . The present study is confined to the pH 
range 8.3 -10.6 and 8.5 - 1 1 for Ca(II) and Mg(II) complexes of Gin, respectively. In 
this pH range the ligand is expected to be present in the form of LH and L'. Hence the 
plausible metal-ligand species are ML 2 , MLH 2+ , ML 2 H + , which are confirmed by 
MINIQUAD75. The equilibria in these pH ranges can be represented as follows : 

MLH 2+ + LHo ML 2 H + + H f (pH = 8.5-10.0) 

ML 2 H + <=> ML 2 + H + (pH = 9.5-11.0) 

The various forms of sue can be given as LH 2 , LH” and L 2 ” in the pH ranges 3.0- 

6.5, 3. 0-7.0 and 4.5-7.0, respectively. The pH range of present study is 3. 0-9.0 where 
the active forms of ligand are LH' and LH 2 . The species confirmed in the pH range are 
ML 2 ”MLH + , ML 2 H”, MLH + . The formation equilibria can be represented as 

M 2+ + LH 2 <=>MLH + + H + (pH = 3.0 - 5.5) 

MLH + + LH' » ML 2 H' + H + (pH = 5.0 - 7.0) 

ML 2 H' o ML 2 ' + H + (pH = 6.5 - 8.0) 

The variation of species concentration with pH is shown in Fig. 1 for some typical 
systems. In glutamine complexes the concentrations of ML 2 and ML 2 H + are almost 

same. But in succinate complexes ML 2 2 ' is the predominant species compared to 
MLH + and ML 2 H‘. Another important observation is that the concentration of free 
metal ion in the presence of Gin is negligible where as there is still considerable 
amount of free metal in the presence of succinic acid. This indicates that Gin is a 
stronger complexing agent than Sue. 

Biological relevance of present study. The polarity at the active-site cavities in 
enzymes/proteins is less compared to that of bulk water. The presence of acetonitrile 
in aqueous solution considerably decreases the dielectric constant and these solutions 
are expected to mimic the physiological conditions. This equivalent solution dielectric 
constant 3 for protein cavities should be generally applicable when it is possible to 
compare ligand binding to the metal ion in protein and in mixed solvent environments. 
The studies carried out on these systems under the present experimental conditions are 
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useful to understand 1) the role played by the active site cavities in biological 
molecules, 2) the type of complex formed by the metal ion, and 3) the bonding 
behaviour of the protein residue with the metal ion. Hence, the various species 
detected in the present study and their concentrations can be extrapolated to the 
biological systems under similar conditions. 



Fig 1- Species distribution diagrams in 40% acetonitrile-water mixture 
A) Ca(II)-G!n, B) Ca(Il)-Suc, C) Mg(II)-Gin, D) Mg(II)-Suc 
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Abstract 

The charge transfer interactions of ten L-aminoacids namely isoleucine, leucine, phenylalanine, 
tryptophan, valine, asparagines, glycine threonine, histidinemonohydrochloride, lysmemono- 
hydrochlonde with chloranil has been studied in 50% aqueous ethanol by conductometric 
technique employing the models of Gutmann and Keyzer and Randhawa and Lakhani. 
Conductometric titration indicates 1 1 stoichiometry of these complexes The equilibrium 
constant K and other thermodynamic functions have been determined and compared with those 
reported in the literature on the basis of spectroscopic and non-spectroscopic methods and 
models 

(Keywords charge transfer/cond u ctan ce/L-am i noac i ds- chloranil systems) 

Introduction 

Charge transfer complexes are of importance as reaction intermediates 1 and in the 
field of drug-receptor binding mechanism 2 . The stability constants of chloranil with 
various acids involving studies on weak interactions were reported 3 ' 4 . In the present 
communication the charge transfer between ten L-amino acids namely isoleucine, 
leucine, phenylalanine, tryptophan, valine, asparagine, glycine, threonine, 
histidinemonohydrochloride, lysine and chloranil by the conductometric method of 
Gutmann and Keyzer 5 ” 7 and a model reported earlier from our laboratory 8 is reported. 

Materials and Method 

Two sets of experiments were performed, one to determine the composition and 
second to determine formation constants of the complexes. In order to determine the 
stoichiometric ratio of CT complex, equimolar stock solutions (5 x 10" 4 M) of D and A 
in 50% aqueous ethyl alcohol, were mixed in the conductivity cell, keeping the total 
volume constant by varying the volumes of D and A. In the second set of experiments, 


*This work was presented at the International Symposium on Trends in Medicinal 
Chemistry and Biocatalysis held at University of Delhi from 26-29 January ,2000. 
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the stock solutions of donors with concentration 5 x 10' 2 M and acceptor with 
5 x lO^M except in case of tryptophan - chloranil system i.e. SxlO 0 M and 5 x 1 O' 5 M 
were prepared in 50% aqueous ethanol, respectively. In order to determine K the 
concentration of D was varied while keeping the concentration of A fixed. The 
conductance measurements were made at three different temperatures i.e. 293, 303 
and 313 ± 0.01K on a Philipps PR9500 conductivity bridge using a conductivity cell 
pre standardized fitted with platinized platinum electrodes The methodology to 
determine conductivity data and stoichiometry 5 ’ 7 , formation constants 8 and 
thermodynamic functions have been reported in the literature and the data so obtained 
are listed in Tables (1-3). 


Results and Discussion 

The present conductometric results on the composition of amino acid - chloranil 
complexes in 50% aqueous ethanol except lysine monohydrochloride conforms to the 
spectrophotometric studies 9 ' 10 . However, in case of lysinemonohydrochlonde, we 
failed to observe a break at 1:2 (D:A), since one of the amino groups of lysine is 
blocked in lysine monohydrochloride. Further, conductometric data (Table 1) reveal 
that there exists no definite relationship between the values of ct p . a M and Ef as well 
as with pKi of amino acids as expected. This may be attributed to small variation of 
pKi (±0.34) and pK .2 (±0.74) values of amino acids. The absence of any definite 
relationship between K, fn K and pd (Table 2) with a p , E? and pK| (Table 1) may also 
be due to the above mentioned reason. A close look at Table 3 reveals that K values 
not only depend on temperature, the dielectric constant and pH of the medium, the 
nature of the method but also on the type of the model employed. One of the reasons 
for the difference between K values determined by spectral and conductance methods 
may be due to the fact that concentration manifests in spectral while activity in 
conductance method. 

The dipole moment of the complex always depends upon the total polarizability of 
the complex. However, in the calculation of the dipole moment of the complex, the 
depolarizability is ignored on the basis that the polarizability is equal to the sum of the 
individual components. But from the depolarizability data in Table 2, it is evident that 
depolarizability increases over that of the additive values of the individual molecules. 
Thus to obtain the exact value of dipole moment, the depolarizability should be taken 
into account also. Further the inspection of Table 2 (columns, 8-10) shows that p d , the 
depolarization constant grows with the growth of temperature which is reverse to the 
behaviour of K values with temperature as expected. 



Table I- Conductance(a) data for 1 1 CT complexes of L-aminoacid donors with chloranil acceptor in aqueous ethanol (50%) at 293K 
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Hydrochloride 
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The present values of K and other thermodynamic functions have been compared 
with the literature values 11 (Table 4) The sign of literature values of A H which is a 
true measure of strength of interaction and AS do not seem to be a typographical 
mistake, since these values do not fit into Gibbs-Helmholtz equation. In view of this 
we do not trust the literature values of thermodynamic functions determined from 
conductance data by Sahai and Singh 11 The equation used by Sahai and Singh 11 is 
quadratic in K and will give two roots. The authors have not given any justification for 
choosing value of K for further calculations and investigations whereas in the model 
employed in the present studies, only one value of K is obtained and hence there is no 
such ambiguity as above for further investigations. Further from the equilibrium 
constant data it is clear that the interaction is stronger in case of trytophan than in case 
of phenylalanine. Theoretical calculation of Pullman reveals that trytophan having an 
indole nucleus should behave as a better electron donor than phenylalanine 12 . 

Table 4- A comparison of thermodynamic functions of tryptophan and phenylalanine * chloraml 


systems 


Donor 

Literature values(a) (295 K) 


Present studies (293 K) 


K 

AH 

AS 

K 

-AH 

AS 

Tryptophan 

1 77 77 

22 99 

79 33 

295 99 

32.00 

-62 01 



(-22 99) 

(-34 88) 




Phenylalanine 

133 33 

22 23 

72 77 

243 68 

11 89 

+5 05 



(-22 23) 

(-34 71) 





a - ref (11). The quantities within parenthesis are the redetermined values 
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Abstract 

Rapid, simple and sensitive spectrophotometric methods for the determination of pyrocatechol, 
dopamine hydrochloride, levo dopa, methyl dopa and adrenaline hydrochloride in either pure 
form or its pharmaceutical preparations is described The first method is based on the 
interaction of catecholamine derivatives with iron(III) and subsequent reaction with 
ferricyamde in presence of hydrochloric acid medium to yield prussian blue coloured complex 
with 7. max of 730 nm In the second method, the interaction of nitrite ions with the 
catecholamines in neutral medium m presence of aluminium ions to yield a red product in 
alkaline medium with I max of 500-510 nm The optical characteristics, interference studies 
and application to pharmaceutical preparations have been reported 

(Keywords catecholamines/spectrophotometiy/alumimum lons/ferricyamde/pharmaceuticals) 

Introduction 

Catecholamine derivatives are aromatic vic-diols in which either the 3- or 4- 
position is unsubstituted and these positions are not sterically blocked. These 
derivatives are used as drugs and have wide application 1 . There are various 
instrumental techniques for the determination of these catecholamine derivatives We 
have already reported different spectrophotometric procedure 2 * 6 for the determination 
of catecholamines. In continuation of our work in this field, we are reporting two 
more sensitive and rapid methods for the spectrophotometric determination of 
pyrocatechol (PCL), dopamine hydrochloride (DPH), levo dopa (LDP), methyl dopa 
(MDP) and adrenaline hydrochloride (ADH) in pure form as well as in its 
pharmaceutical formulations. Optical characteristics like Beer's law range, molar 
absorptivity, Sandell's sensitivity etc , have been reported. Commonly used excipients 
do not interfere in the present methods. Reaction mechanisms are discussed. 
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Materials and Method 

A JASCO model UVIEC-610 spectrophotometer with 1.0 cm matched cells was 
used for absorbance measurements. Pharmaceutical grade PCL, DPH, LDP, MDP and 
ADH were purchased from Sigma Chemical Co. (USA). Aluminium chloride was 
purchased from Merck (Darmstadt, Germany). All other chemicals and solvents used 
were of analytical reagent grade. Commercial dosage forms were purchased from 
local sources. Deionized water was used throughout. 

Stock solutions of catecholamine derivatives were prepared by dissolving 0.05 g 
of each drug in 100 ml of water. The working standard solutions of the drug were 
prepared by further dilution. A 1 % aqueous sodium nitrite, 0.2% aluminium chloride, 
1M sodium hydroxide, 0.002M potassium ferricyanide, 0.1 M ferric chloride and 5M 
hydrochloric acid solutions were used for the experiments. 

General procedure : Method A : Aliquots of standard solutions of PCL (2.5-140 
Mg), DPH (6.25-275 Mg), LDP (2.5-100 Mg), MDP (6.25-500 Mg) or ADH (5-650 Mg) 
were transferred into a series of 25 ml calibrated flasks. To each of the flasks, 1.0 ml 
of 1% sodium nitrite and 0.5 ml of 0.2% aluminium chloride were added. After 5 min. 
2.0 ml of 1M sodium hydroxide was added, made upto the mark and shaken well. The 
absorbances of the red coloured products were measured at 500-5 1 0 nm against a 
reagent blank and calibration graphs were constructed. Optimum ranges are given in 
Table 1. Method B : Aliquots of standard solutions of PCL (2.5-80 Mg), DPH (2.5-80 
M§)> LDP (2.5-60 Mg), MDP (1.25-25 Mg) and ADH (2.5-80 Mg) were transferred into 
a series of 25 ml calibrated flasks. Requisite amounts of potassium ferricyanide and 
iron(III) were added as shown in Table 1 and left for 5 min. Then, 1.0 ml of 5M 
hydrochloric acid was added and diluted with water upto the mark. The contents were 
mixed thoroughly and the absorbance values were recorded at 730 nm against the 
reagent blank. Calibration graphs were constructed. 


Table 1- Experimental conditions. 


Catechol 

derivatives 

1 % NaN0 2 

Optimum reagent concentrations 

0 2% A1C1 3 IMNaOH Iron(Ill) 

Potassium 

ferricyanide 

PCL 

*1 ml *(0.5-3 
ml) 

0.5 ml 
(0.3-3ml) 

2mi (i-4ml) 

0.05 M 

5ml (1-7 ml) 

0.001 M 4 ml 
(1-5 mi) 

DPH 

1 ml (0.5-3 
ml) 

0.5 ml 
(0 3-3ml) 

2ml(l-4ml) 

O.i M 

3 ml (2-8ml) 

0 002 M 3 ml 
(2-8 mi) 


Table 1 Contd 
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LDP 

I ml (0 5-3 

0 5 ml 

2m 1 (l-4ml) 

0 05 M 

0.01 M 


ml) 

(0 3-3ml) 


4m! (1-5 ml) 

4 ml (1-5 ml) 

MDP 

1 ml (0 5-3 

0 5 ml (0 3- 

2ml (l-4mi) 

0 ! M 

0 002 M 


ml) 

3m!) 


5 ml (1 -7ml) 

5 ml (1-7 ml) 

ADH 

1 ml (0 5-3 

0 5 ml 

2ml(l-4mi) 

0 05 M 

0 001 M 


ml) 

(0 3-3ml) 


4 ml (!-5rn!) 

4 ml (1-5 ml) 


“Used in the proposed procedure 

b Range for maximum absorption and stability (values in parenthesis) 


Procedure for the assay of pharmaceutic samples Fifteen tablets were weighed 
and finely powdered. The powder amount equivalent to 50 mg of LDP or MDP (for 
injection of PH or AH, an appropriate volume of the sample) was dissolved in water 
and filtered. The fi iterate was made upto 1 00 ml and an appropriate aliquot of the drug 
solutions were treated as described in the general procedure. 

Results and Discussion 

The optical characteristics and precision data for both the methods are given in 
Tables (2 & 3). Method A involves the reaction of catechol derivatives with sodium 
nitrite in neutral medium in presence of aluminium ions to produce an intense yellow 
coloured product which turns to deep red after addition of NaOH. It is presumed that 
anhydrous aluminium chloride not only acts as lewis acid for the nitration of 
catechols, but in higher concentration, it stabilizes the nitro derivatives by complex 
formation, which readily gets decomposed in alkaline medium. A bathochromic shift 
from 420 nrn to 500-5 1 0 nm in alkaline medium is noticed. The probable mechanism 
is shown in Scheme 1. Method B involves the reduction of iron(III) in aqueous 
medium to form iron(ll), which subsequently chelate with ferricyanide to form 
Prussian blue coloured product. 

Both the reactions .are completed within 5 min. at room temperature (26±1°C). 
The proposed methods were found to be free from the interference of commonly 
encountered excipients and additives such as sodium chloride, talc, starch, stearic 
acid, gelatin, gumacacia, glucose, lactose, sucrose and sodium alginate. Vitamin-B6 
and vitamin-C (upto 4000 ppm) do not interfere in method A, while any reducing 
agent interferes in method B. Some typical analysis of the results of tablets and 
injections are presented in Table 4. 
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Table 2- Optical characteristics and precision data tor method A 


Characteristics / Parameter: 

PCL 

DPH 

LDP 

MDP 

ADH 

Colour 

Red 

Red 

Red 

Red 

Red 

'-max (nm) 

510 

505 

500 

500 

500 

Stability (h) 

20 

22 

18 

10 

10 

B-L range (jig/ml) 

0 1-6 4 

0 25-1 1 

0 1-6 8 

0 25-20 

0 2-26 

LOD (jig/mS) 

0 0544 

00918 

0 0379 

02148 

0 1523 

LOQ (pg/rn I) 

0 1815 

0 3059 

0 1266 

0.7161 

0 5077 

Optimum photometric raige 
( jig/ml) 

0 5-5 2 

0 75-10 

0 5-5 5 

0 75-18 

0.8-24 

Molar absorptivity (1 mol' 1 cm' 1 ) 

1 35x 1 0 4 

9 4\10 3 

2 1 x 1 0 4 

6 45x1 0 3 

6 2xl0 3 

Sandell's sensitivity (pg/cm : ) 

0 0085 

0 0162 

0 0094 

0 032 

0 029 

Regression equation (Y)* 

Slope (a) 

0 04798 

0 0568 

0 1031 

0 0304 

0 0257 

Intercept (b) 

-0 0075 

0 0049 

-0 0105 

-0 0013 

00138 

Correlation coefficient (r) 

0 9986 

0 9987 

0 9969 

0 9998 

0 9995 

RSD % (n=5) 

0 84 

1 1 

1 2 

1 3 

1 4 

% Range of error (n=5) at 95% 
confidence limit 

±0 30 

±0 42 

±0 38 

±0 50 

±0 52 


*Y= ax+b, where jc is the concentration of catecholamines in jig/ml 

r -\ 



Nitroso derivative 



Nitro catechols (Red) 

(yellow) 

where R = H-+PCL . -CH 2 CH 2 NH 2 (DPH), -CH 2 C(CH 3 )-NH 2 COOH->(MDP), 
-CH 2 -CH(NH 2 ).COOH->(LDP) 1 -CH(OH)-CH 2 NHCH 3 ~>(ADH) 

Scheme 1 -Reaction sequence for method A 
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Table 3- Optical characteristics and precision 

data for method B 



Characteristics / Parameters 

PCL 

DPH 

LDP 

MDP 

ADH 

Colour 

Blue 

Blue 

Blue 

Blue 

Blue 

uw 

730 

730 

730 

730 

730 

Stability (h) 

1 

1 

1 

1 5 

1 

B-L range (pg/ml) 

0 1-3 2 

0 1-2 0 

0 1-2 4 

0 05-1 0 

0 1-3 2 

LOD (pg/ml) 

0 0501 

0 0335 

0 0301 

0 0026 

0 0456 

LOQ (pg/ml) 

0 1671 

0 1116 

0 1004 

0 0087 

0 1523 

Molar absorptivity (1 mol' 1 cm' 1 ) 

5 3xl0 4 

7 2xl0 4 

1 45xl0 5 

1 14x10 s 

4 4xl0 4 

Sandell's sensitivity (pg/cm 2 ) 

0 0020 

0 0026 

0 0014 

0 0018 

0 0042 

Optimum photometric range 
(pg/ml) 

021-28 

0 2-1 8 

0 21-2.0 

0.1-0 8 

0 2-2 8 

Regression equation (Y)* 






Slope (a) 

0 2087 

0 2733 

0 2604 

0 5013 

0 1716 

Intercept (b) 

0 0314 

0 0241 

0 0772. 

-0 0053 

0 0238 

Correlation coefficient (r) 

0 9971 

0 9796 

0 9836 

0 9980 

0 9959 

RSD % (n=5) 

0 5746 

0 747 

0 7699 

0.2795 

0 3462 

% Range of error (n=5) at 95% 
confidence limit 

±0 7977 

±1 038 

±1 068 

±0 3882 

±0 4805 


*Y = ax + b, where * is the concentration of catecholamines in pg/ml 

Table 4- Determination of catechol derivatives in pharmaceutical preparations 


Drug Label claim % Recovery ± % RS D * 

BP Method 7 Method A Method B 

Injections 

DPH* 200/5 ml 

DPH b 200/5 ml 

_ ~ ‘ ~ Table 4 Contd 


98 85±1 25 99 80 ± 0 40 99 90 ± 0 42 

99 20±0 82 99 75 ± 0 50 99.85 ± 0.75 
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Table 4 Contd 


ADH C 

200/5 ml 

99 40 ± 0 55 

98 95 ± 0 50 

98 85 ± 0 65 

ADH d 

200/5 ml 

100 00 + 0 90 

99 6 ± 0 80 

99 4 ± 0 95 

Tablets 





LDP C 

500 

9878 ± 1 10 

99 2 ± 0 70 

99 80 ± 0 90 

MDP f 

250 

98 78 ± 1 10 

98 30 ± 1 3 

98 84 ± 1 02 


*n =5 , Marketed by • *TTK Pharma, b Troika Paren; 'Harson Laboratories, d Wallace Pharma, e Sun 
Pharma; f Mermd Pharma 


Conclusions 

The proposed methods are found to be simple, rapid and sensitive than most of the 
available methods in literature. Heating or extraction is avoided and compounds like 
resorcinol, phloroglucinol, hydroquinone and pyrogallol do not give the colour 
reaction for method A. The statistical parameters and the recovery study data clearly 
indicate the reproducibility and accuracy of the methods. Thus the methods can be 
adopted as an alternative to the existing methods. 
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Abstract 

A simple, convenient and highly sensitive method of determination of micro amount of nickel 
has been developed The coloured complex formed on reaction of di(o-ethylphenyl) carbazone 
with nickel(II) to give 1 . 2 complex in presence of 1,10-phenanthrohne It is found that the 
extractability of the nickel di(o-ethylphenyl) carbazonate into chloroform can be increased by 
the addition of phenanthrohne Synergistic effect for nickel di(o-ethylphenyl) carbazonate with 
phenanathroline has been investigated. 

(Keywords • spectrophotometry/synergistic effect/mckel adducts) 

Introduction 

In continuation of our studies on substituted diphenylcarbazones and dithizone as 
analytical reagent 1 * 3 , we have examined the synergic extraction of nickel(II) chelate 
employing bidentate chelating agent di( o-ethylphenyi) carbazone in the presence of 
an adduct forming substance, 1,10-phenanthroline has received the considerable 
attention for the spectrophotometric determination of nickel(II) 4 ' 5 . The method is 
elegant, sensitive and simple for extractive-photometric determination of micro 
amount of nickel. 

In our present investigation it has been found that the extractability of the 
nickel(II) di( o-ethylphenyl) carbazone complex, Ni(OEPC) 2 , into chloroform can be 
increased by the addition of phenanthroline and the extracted complex is stable for 6 h 
which is sufficient to perform analysis. Such synergistic effects are due to an increase 
in the distribution ratio of a metal chelate on addition of a base such as phenanthroline 
to the metal chelate. The detailed mechanism of synergistic effects for various metal 
chelates in presence of neutral bases has been investigated extensively 6 and a few 
applications to analytical chemistry have also been reported 7 . In general, the 
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mechanism seems to involve the formation of a base adduct of a metal chelate. The 
synergistic effect in the present extraction results from the formation and the 
preferential extraction of a phenanthroline adduct of Ni(OEPC) 2 . 

Materials and Method 

The reagent di(o-ethylphenyl) carbazone was prepared by the reported method 8 
involving persulphate oxidation of the corresponding carbazide and then purified 9 . 
Chloroform (B.D.H.) was purified by the method suggested by Vogel 10 . The constant 
boiling fraction was collected and used. Nickel perchlorate and 1 , 1 Q-phenanthroline 
(Fisher Reagent Grade) were used as received. The absorbances were measured by 
using the Baush and Lomb Spectronic 2000 recording spectrophotometer. 

To 10ml of a solution containing from 0.! to lOpg of nickel (upto 25 pg may be 
accommodated) were added 5m! of phthalate or acetate buffer of pH 6.0 (or dilute 
ammonia may be used) followed by 15ml of a chloroform solution containing 7xlQ‘ 5 M 
Ni(OEPC) 2 and 3xlO' 5 M phenanthroline. The two phases were shaken for 5 min. 
phases were separated after allowing them to settle. The chloroform layer was back 
extracted with 10ml of 0.1 M sodium hydroxide shaking vigorously for about 1 min. 
The organic layer was separated and the absorbance of the extract was measured at 
515nm against the reagent blank. The nickel content in unknown solutions was 
computed from the standard calibration curve. 

Typical results are tabulated below : 


Ni taken jig 

Absorbance at 515nm 

Ni found, jig 

0.71 

0.13 

0 73 

1.46 

0.27 

1.47 

2.28 

0.42 

2.31 

3.04 

0.56 

3.01 

3.37 

0.62 

3.29 

3.86 

0.71 

3.87 

4.84 

0.89 

4.91 
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Results and Discussion 

The nickel(II) di( o-ethylphenyl) carbazonephenanthroline adduct system exhibits 
maximum absorbance at 515nm when measured against the reagent blank and obeys 
Beer's Law over the rage of 0.1-0.25mg of Ni". The Sandell's sensitivity of the colour 
reaction is 0.00049|ig cm" 2 and molar absorptivity 5.33xl0 4 lit. mole' 1 cm' 1 at 515nm. 

The results of extraction studies of nickel(II) di(o-ethylphenyl) carbazone system 
shows that nickel can be quantitatively extracted into chloroform in the pH range 6.0- 
11.0. Below and above the pH range, the extraction of nickel was incomplete. The 
extracted species exhibited constant and maximum absorption, when the pH of the 
aqueous phase was maintained with this range. 

Various solvents such as benzene, toluene, amyl alcohol, carbon tetrachloride and 
chloroform were tested for the extraction of nickel complex, but chloroform was 
found to be the most effective. With 10ml of chloroform, the complex was recovered 
to the extent of 99.98% by a single extraction, whereas with other solvents, it was 
found to be less than 95%. For complete extraction of the complex, 5 min. shaking 
was adequate. 

The mixed ligand complex is sufficiently stable to permit the removal of excess of 
di(o-ethylphenyl) carbazone by back - extraction with 0.1M sodium hydroxide, so a 
'monocolour’ method is applicable. The molar absorptivity of this complex is 5.33xl0 4 
I mole^cm' 1 which makes a method based on this reaction about four times as 
sensitive as the dimethylglyoxime method used as oxidising agent 11 . 
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Abstract 

Three simple, rapid and cost-effective methods based on titnmetry and spectrophotometry' are 
described for the determination of albendazole in pure form and in dosage forms In titnmetry, 
the drug is dissolved in glacial acetic acid and titrated with acetous perchloric acid with visual 
and potentiometnc end point detection, crystal violet being used as indicator for visual titration 
Spectrophotometry involves adding different amounts of drug to a fixed amount of perchloric 
acid-crystal violet mixture and measuring the increase in absorbance at 570 nm Both the 
titrimetric methods are applicable over 2-20 mg range, and m spectrophotometry, Beer's law is 
obeyed over range 5-60 pg ml* 1 The apparent molar absorptivity is calculated to be 2 18 x I0 3 1 
mol' 1 cm' 1 with a Sandell sensitivity of 121 9 ng cm' 2 The limits of detection and quantification 
are calculated to be 1.02 and 3.39 pg ml' 1 , respectively The methods when applied to the assay 
of albendazole m tablet gave satisfactory results (97 26 to 102 53 %). The accuracy and 
reliability of the methods were further ascertained by parallel determination by a reference 
method and by recovery studies via standard-addition technique 

Keywords : Albendazole/determination/titrimetry/spectrophotometry/non-aqueous medium/ 
formulations 


Introduction 

Albendazole, chemically, 5-propyl thio-lH benzimidazole methyl carbamate 1 (Fig. 
1) which in therapeutics is known as anthelmentic drug having a wide spectrum of 
activity 2 . The therapeutic importance of albendazole justifies research to develop 
analytical methods for its determination in biological samples and in pharmaceuticals. 
The most widely used technique for the assay of albendazole has been the high 
performance liquid chromatography (HPLC) and has been used for the determination 
of the drug in biological fluids 3,4 , plasma 5 ' 9 , veterinary formulations 10 and meat 
samples 11 . Ion-pair liquid chromatography with uv-detection has been applied to the 
determination of the drug residus in milk and drug metabolites 12 , and in cattle feed 13 . 
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The same technique with fluorescence detection has been used for drug in milk 14 . 
Liquid chromatography (LC) with fluorescence detection 15 and uv-detection 16 has 
found applications in the quantification of albendazole residues in cattle liver and in 
pork muscle tissue, respectively. Even-hyphenated techniques such as LC-MS and 
GC-MS have found use in the determination of drug residues in food products 17 and in 
cattle liver 18 , respectively. 


H 



Fig 1- Structure of albendazole 

Not many methods are found in the literature for the determination of albendazole 
in pharmaceutical products. Five methods based on the use of HPLC 19 ' 23 have been 
proposed for bulk drug and for formulations and are less sensitive with linear ranges 
being 20-100 pg ml' 1 19 0. 1-0.4 mg ml' 1 20 and 50-300 pg ml' 1 2! . The drug in tablets 
has been assayed by voltammetry 24 and uv-spectrophotometry 25,26 . Two visible 
spectrophotometric methods, one based on ion-pair complex formation between the 
drug and acid dyes such as bromocresol green, bromophenol blue, bromophenol red 
and bromothymol blue followed by extraction into chloroform and absorbance 
measurement at 420 nm has been proposed by Sane et. al. 21 A similar reaction but 
involving picric acid 28 as the ion-pair reagent has served as the basis for the 
spectrophotometric determination of albendazole. Zarapkar and Deshpande have 
reported a redox method using Folin-Ciolteatu reagent 29 for the assay of albendazole 
in tablets and syrup. Another redox method involves the treatment of the drug with a 
known excess of N-bromosuccinimide (NBS) 30 followed by the determination of the 
surplus oxidant spectrophotometrically based on chrge-transfer complex reaction or 
using the dye Celestine Blue One of the present authors 3 ' has recently reported three 
procedures using bromate-bromide mixture and methyl orange dye, based on 
oxidation-bromination reaction, for the assay of drug both in pure form and in 
formulations. Though the procedures are fairly sensitive, they suffer from such 
disadvantages as extraction step 27,28 , lack of selectivity 29 ’ 31 , and use of an unstable 
reagent’ 0 . 

The present paper describes three methods based on the basic property of the drug 
molecule. In titrimetric procedures, the drug solution in glacial acetic acid is titrated 
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directly with acetous perchloric acid, the end point being determined either visually 
using crystal violet indicator or potentiometrically using combined glass-modified 
SCE system Spectrophotometry involves the treatment of a fixed amount of 
perchloric acid-crystal violet mixture with drug and measurement of absorbance at 
570 nm. The increase in absorbance is related to drug concentration. 

Materials and Method 

Apparatus All absorbance measurements were made with a Systronics Model 106 
digital spectrophotometer provided with matched 3 -cm quartz cells. Potentiometric 
titration was performed with an Elico-120 digital pH meter provided with a combined 
glass-SCE system, the KCI of the salt bridge being replaced with 0.1M lithium 
perchlorate solution in glacial acetic acid. 

Reagents and standards All chemicals used were of analytical reagent grade. All 
solutions were made in glacial acetic acid unless specified otherwise. 

Perchloric acid (0 01 M) To 4.5 ml of 70 % perchloric acid (S. d. Fine Chem, 
Mumbai, India) was added 150 ml of glacial acetic acid, mixed well; added 10.5 ml of 
acetic anhydride and allowed the solution to cool for 30 min; finally diluted to 500 ml 
with glacial acetic acid and allowed to stand overnight. This perchloric acid (~0.1 M) 
was diluted to get 0.01 M acid with glacial acetic acid and standardized with pure 
potassium biphthlate and crystal violet indicator. 

Crystal violet indicator. Prepared by dissolving 0.1 g of dye in 100 ml glacial 
acetic acid. 

Perchloric acid-crystal violet mixture (1.5 mM HC10 4 - 0.25 mM crystal violet) 
Prepared by mixing 15 ml of 0.01 M perchloric acid and 10 ml of 1000 pg ml' 1 crystal 
violet solutions and diluting to 100 ml with glacial acetic acid. 

Standard drug solution Pharmaceutical grade albendazole was procured from 
Cipla India Ltd., Mumbai, as gift, and was used as received. A stock standard solution 
containing 2 mg ml' 1 albendazole was prepared by dissolving 500 mg of pure drug in 
glacial acetic acid and diluting to the mark in a 250 ml calibrated flask. This solution 
(2 mg ml' 1 ) was used for titrimetric work, and for spectrophotometric work, the same 
was diluted appropriately with glacial acetic acid to get 100 pg ml' 1 working 
concentration. 

Procedures 


Visual titration (Method A). A 10 ml aliquot of standard drug solution containing 
2-20 mg of albendazole was pipetted out into a clean and dry 100 ml titration flask, 2 



274 


K BASAVAIAH et al 


drops of crystal violet indicator was added and titrated with standard 0 01 M 
perchloric acid to an emerald green end point. The amount of drug in the measured 
aliquot was calculated from 

Amount (mg) = VMR 

where V = volume of perchloric acid required, ml 
M = relative molecular mass of drug, 

R = molarity of perchloric acid. 

Potentiometric titration (Method B) A 10 ml aliquot of standard drug solution 
equivalent to 2-20 mg of albendazole was pipetted out into a clean and dry 100 ml 
beaker and the solution was diluted to 30 ml by adding glacial acetic acid. A 
combined glass-SCE (modified) system was dipped in the solution. The contents were 
stirred magnetically and the titrant (0.0 1M HCI0 4 ) was added from a micro burette. 
Near the equivalence point, the titrant was added in 0.2 ml increments. After each 
addition of titrant, the solution was stirred magnetically for 30 S and the steady 
potential was noted. The addition of titrant was continued until there was no 
significant change in potential on further addition of titrant. The equivalence point 
was determined by applying the graphical Gran's plot method. The amount of drug in 
the measured aliquot was calculated as described under visual titration. 

Spectrophotometric assay (Method C) Different aliquots (0 5 - 6.0 ml) of 
standard 100 pg ml' 1 drug solution were accurately transferred into a series of 10 ml 
calibrated flasks. An exactly measured (2 ml) perchloric acid-crystal violet mixture 
was added to each flask and the volume was diluted to the mark with glacial acetic 
acid, mixed well and absorbance was measured at 570 nm against a reagent blank. 
The increasing absorbance values at 570 nm were plotted against the concentration of 
drug to obtain the calibration graph. The concentration of the unknown was read from 
the calibration graph or calculated from the regression equation obtained from Beer's 
law data. 

Procedure for tablets. Twenty tablets were weighed and ground into a fine 
powder. An amount of powder equivalent to 500 mg of albendazole was weighed 
accurately into 250 ml calibrated flask, 150 ml of glacial acetic acid added and shaken 
for about 20 min. Then, the volume was made up to the mark, mixed well and filtered 
using a quantitative filter paper. First 10 ml portion of the filtrate was discarded and a 
suitable aliquot was then subjected to analysis by titrimetry The filtrate equivalent to 
2 mg ml 1 was diluted appropriately to obtain 100 pg ml" 1 drug solution and analysed 
by spectrophotometry as described under the general procedure. 
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Results and Discussion 

Several titrimetric 31 , voltammetric 24 and spectrophotometric 30,31 methods have 
been described for the determination of albendazole based on latter's redox property. 
A couple of extraction spectrophotometric methods 27,28 utilizing the tertiary nitrogen 
of the drug molecule have also been proposed. The present methods are based on the 
basic property of albendazole employing two techniques. The methods are based on 
the principle that substances, which are too weakly basic in aqueous medium, exhibit 
enhanced basicity in non-aqueous solvents thus allowing their easy determination. In 
the present titrimetric methods, the weakly basic property of albendazole was 
enhanced due to the non leveling effect of glacial acetic acid and titrated with 
perchloric acid with visual and potentiometric end point detection. Crystal violet gave 
highly satisfactory end point for the concentrations of analyte and titrant employed. A 
steep rise the potential was observed at the equivalence point with potentiometric end 
point detection (Fig 2). The Gran's plot method was applied to ascertain the 
equivalence point (Fig 3). With both methods of equivalence point detection, a 
reaction stoichiorretry of 1: 1 (drug: titrant) was obtained which served as the basis for 
quantitation. Using 0.01 M perchloric acid, 2-20 mg of albendazole was conveniently 
determined. The relationship between the drug amount and the titration end point was 
examined. The linearity between two parameters is apparent from the correlation 
coefficients of 0.9983 and 0.9996 obtained by the method of least squares for visual 
and potentiometric methods, respectively. From this, it is implied that the reaction 
between albendazole and perchloric acid proceeds stoichiometrically in the ratio 1 : 1 
in the range studied. 



Fig. 2- Potentionetric end "point point. 
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Crystal violet is a dye exhibiting violet colour in the base form and emerald green 
in the acid form. The present spectrophotometric method is based on the facts that the 
colour of the dye is dependent on the pH of the solution and that the colour change is 
not sudden but occurs continuously as the pH changes over a definite range. To a 
fixed amount of acid-dye mixture where the dye is in the acid form (emerald green), 
different amounts albendazole were added. This caused a progressive increase in pH 
of the solution because of neutralization of acid by the added drug (base), and as a 
result, the concentration of the acid form of the dye increases. This is shown by the 
proportional increase in the absorbance of the solution at 570 nm (Fig 4) which is 
corroborated by the correlation coefficient of 0.9982. 



Fig. 4-Calibration graph 
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In a preliminary study, 20 pig ml' 1 crystal violet in basic medium was found to 
exhibit a convenient absorbance at 570 nm. In the presence of 2 ml of 1.5 mM 
perchloric acid in a total volume of 10 ml, this absorbance decreased to a constant 
minimum. Hence, different amounts of drug were treated with a fixed amount of acid- 
dye mixture i.e, 2 ml of 1.5 mM HC!0 4 - 0.25 mM crystal violet (100 pig ml' 1 ) to 
determine the concentration range of the drug that could be determined by the method 
of absorbance transitions of the dye accompanying the pH changes. The dye colour 
was found to be stable for several hours, and the order of addition of reactants was not 
critical. 

The increasing absorbance values at 570 nm were plotted against the increasing 
concentration of drug to obtain a calibration graph. Beer's law is obeyed over the 
concentration range 5-60 jag ml' 1 , the equation of the line being 

A= 0.013 + 0.007 C 

where A is absorbance and C concentration in pig ml' 1 . The correlation coefficient of 
the calibration plot was calculated to be 0.9982 (n=7) confirming a linear increase in 
absorbance with increasing concentration of albendazole. The calculated molar 
absorptivity was found to be 2.2 x 1 0 3 1 mol' 1 ml' 1 at 570 nm and the Sandell 
sensitivity was 122 ng cm' 2 . The limits of detection and quantification were calculated 
form the standard deviation of the absorbance measurements obtained from a series of 
seven blank solutions. The limits of detection and quantification established according 
to IUPAC definitions 32 were 1 .02 and 3.39 pig ml' 1 , respectively. 

Accuracy and precision 

The accuracy and precision of the methods were established by analyzing the pure 
drug solution at three different levels. The relative error (%), which is a measure of 
accuracy, is < 2% revealing high accuracy of the methods. The coefficient of 
variation, which is an indicator of precision, is less than 1.5 % and speaks of excellent 
precision of the methods. The results of the study are complied in Table 1. 

Application 

The proposed methods were successfully applied to determine albendazole in 
tablets. The same batch tablets were analysed by an established procedure for 
comparison. The results obtained by the proposed methods agree well with these of 
reference method and with the label claim. 



Table 1- Evaluation of accuracy and precision of the methods 
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The results were also compared statistically by a Student's t-test for accuracy and 
a variance ratio F-test for precision with those of the reference method 33 at 95% 
confidence level as summarized in Table 2. The results showed that the calculated t- 
and F- values did not exceed the tabulated values inferring that proposed methods are 
as accurate and precise as the reference method 33 . 


Table 2- Results of assay of tablets by the proposed methods 


Tablet brand 

name* 

Nominal value 

(mg) 

Method A'* 1 

Method B $ 

Method C ,; 

Official 

method 

Alminth a 

200 

99 86 ± 1 25 

100 24 ±0 68 

102 36 ± 1 54 

100.98 ±0 92 



t=l 63 

t= 1 26 

t=l 78 




F=1 88 

F=1 79 

F=2.86 


Albazole b 

400 

100 76 ± 1 06 

100 42 ± 1 18 

97 63 ± 0.85 

99 32 ±0.70 



t=2 58 

t=l 60 

t=1.41 




F=2 29 

F=2 84 

F-l 47 


Dispel 0 

400 

98 62 ± 1 36 

100 63 ±0 58 

99.28 ± 1 16 

100 28 ±0 91 



t~2 31 

t=0 64 

1=1.52 




F=2 23 

F=2.46 

F=1 62 


Nopar d 

50 

99 83 ± 0.97 

99 88 ± 1.28 

102.28 ± 0.82 

101.36 ±1 84 



t=l 72 

t=l 30 

t=I 09 




F-3.59 

F-2.07 

F=5.03 



♦Marketed by a. Torrent Pharmaceuticals, 
b • Wings Pharma; c Indchemie; d • Malladi Drugs , 

Tabulated t-value at 95% confidence level is 2.77 A & C, 2 37 for method B 
Tabulated F-value at 95% confidence level is 6.39 A & C, 9.28 for method B. 
v|i ‘ Mean value of five determinations 
$ : Mean value of three determinations 



Table 3- Results of recovery study by standard addition method. 
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Accuracy and reliability of the methods were further ascertained by performing 
recovery experiments. To a fixed amount of the drug in tablet (pre-analysed), pure 
drug at three different levels was added, and the total was found by the proposed 
methods. Each test was repeated three times. The results compiled in Table 3 show 
that the recoveries were in the range 96.8 - 103.2% indicating that commonly added 
excipients to tablets such as talc, starch, gelatin, sodium alginate, magnesium stearate, 
calcium gluconate, calcium dihydrogenorthophosphate etc, did not interfere in the 
determination. 

Conclusions 

Albendazole has been determined in tablets using three different techniques. The 
developed methods are rapid, simple and accurate. The proposed spectrophotometric 
method is comparable in sensitivity to many of the existing methods including HPLC 
procedures. The procedure is free from tedious steps like extraction or heating and 
involves least number of experimental variables which is reflected in high precision. 
All the three methods are applicable over long dynamic concentration range and can 
serve as useful reference methods which could be used for a specific albendazole 
assay. 
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Abstract 

In the present paper an attempt has been made to obtain and discuss the growth and decay of 
weak discontinuities m characteristic plane for reacting gases, when chemical reaction between 
rt species are taken as the source of non-equilibnum effects The critical amplitutde of the 
initial compressive-disturbance has been determined and it is shown that any compressive- 
wave with an initial amplitude greater than critical one will terminate into a shock-wave, while 
compressive-wave with an initial amplitude less than the critical one will result m a decay 
With the help of Laplace-transformation exact solutions and solution upto first order of small 
quantities have been obtained and discussed for expansion-wave 

(Keywords weak-wave/non-equlibnum flow/reactmg gases) 

Introduction 

Bauer and Bass 1 have demonstrated that gas maintained in ambient vibrational 
and, radiative non-equilxbnum by thermal radiation can amplify acoustic-waves 
Using gas dynamic approach Shrinivasan and Vincenti 2 have obtained a criteria for 
acoustic instability in a gas with ambient vibrational and radiative-non-equijibrium. 
Ram and Smgh 3 have discussed the growth and decay of weak discontinuities in 
relativistic fluids with vibrational-relaxation Johnnesen 4 has analysed the vibrational 
relaxation region by Rayleigh method. Ram and Pandey 5 have discussed the growth 
and decay of acceleration-waves in transient gas flows with vibrational-relaxation. 
Wegener et al. 6 have discussed weak-waves in relaxing flows with a single relaxing 
process. Sharma 7 had obtained explicit critena for the breaking or non-breaking of 
acceleration-waves in plane, axially or spherically symmetric fow of a relaxing gas by 
taking into account a general relaxing gas model. Sharma et al * have discussed the 
behaviour of discontinuity at the wave head propagating through a relaxing gas. Clark 
et al 9 have discussed the dynamics of relaxing gas. Rarity 10 has investigated the 
break-down of characteristic solutions in flow with vibrational-felaxation. Moore and 
Gibson 11 have studied the propagation of weak-disturbances in a gas subject to 
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relaxation effect Shuler 12 has analysed the relaxation processes m multistate systems 
Shaima 13 has investigated the growth and decay properties of first order 
discontinuities headed by wave-front of arbitrary shape m non-equilibrium flows of 
relaxing gas. Borer 14 has discussed the method of characteristics for the computation 
of supersonic flow field m a reacting gas under equilibrium condition Chu 15 has 
discussed the wave propagation in reacting mixture Coleman and Gurtin 16 have 
discussed growth and decay of discontinuities in fluids with internal state variables 
Pandey and Saxena 17 have discussed weak-waves in dissociating gases 

Having such rich back-ground of non-equilibrium flows, in this paper an attempt 
has been made to discuss growth and decay of weak-discontinuity for reacting gases, 
where chemical reaction between n-species has been taken as the source of non- 
equilibnum effect. 

Equation of Motion in Characteristic Plane and Behaviour at Wave-Front 

Neglecting viscosity, heat-conduction and body forces, equations governing one- 
dimensional motion of reacting gases are given by 


U, t +BU, x +D = 0 


( 1 ) 


where U is a column matrix of order (n + 3) x 1 with elements p, u, s, C\, c 2 , . ,c n .B is 

a square matrix of order n+3, with elements b ]2 = pal, b n = —,b„= u and other 

P 

elements zero, D is a column matrix of order (n+3) x 1 with elements 


4 =«ltfe-4 .+ p, 

'=> V 


, vwp , 

X J 


^,=—0 = 1 , .,«) 

1 t=l P 


^ 2=0 


P> K P> s , T being density, velocity, pressure entropy and temperature of gas, ctf is 
frozen sound speed and A h c 2 and co ; are affinity, concentration and rate of production 
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of Uh species v = 0, 1,2 for planer, cylindrical and spherical symmetry A comma 
followed by an index denotes partial denvatives with respect to that index 0 is a null 
matrix of order n + 3 and c, is material derivative (i taking values 1 to «). 

A function U(x,t) satisfying equation (1), every where except at the characteristic 
curve YAf) is said to be a weak-discontmuity at this curve and propagate along the 
characteristics, where U is continuous and U, t and U, x may suffer a finite jump. There 
are («+ 3) families of characteristics of equation two of which 


dx 

— = u±a, 
dt 1 


( 2 ) 


represent the wave propagation mix direction with an effective local-speed of sound 
a/ and 


dx 

— = u 
dt 


(3) 


repeated for (n+ 1) times represent the particle path. 

In studying the wave phenomenon governed by hyperbolic equations, it is usually 
more natural and convenient to use the characteristics of governing system as the 
reference coordinate system Following Chu 15 and introducing two characteristic 
variables a and vj/, equations of motion in characteristic plane are given by 

(p, a - a f p u, a .) t,* - p,y t, a = 0 (4) 


fj> \ 


Ts ,o.~ U,« = 0 

(5) 

V*=l J 


C Ka-C,t a =0,(l= 

(6) 


P,v + p d f u^+ aU £ -zr4+P Ci 


i=i' 






(V 
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The boundary conditions at the characteristic wave front a = 0 are 

P.v =« v = P,v = c,, v = 0.(/=l,. 

= 1 

P,a — Po a fo u ,a~ O' = ^ ) 

Xia = 0 

*, v a f0 

where flow variables associated with the undisturbed medium ahead of the wave 
signified by the subscript 'o'. 

At wave-front, a = 0 

u = ^ — def a. 

^/o’0 a 


Introducing dimensionless parameters 


^ = (T 0 + \)a * vj/ * 

fundamental equation for the evaluation of a weak-wave reduces to 
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dZ 

dx\ 


a + 


2r\ + lJ 


Z+X5 2 =0 


Solution of equation (8) is 

5 = [(l + 2T 1 ) W2 e ar '{l + XIF(T 1 )}]- 1 , 


( 8 ) 


( 9 ) 


where 



g-(an) 

(1 + 2ti) v/2 


■dr\ 


and 


r=i +(«?)„ 

The function W(r |) plays an important role in the breakdown of the characteristic 
solution. As W(t\) > 0, only compressive-wave front may terminate into a shock-wave. 

(i) Plan-Wave-Front 


For plane-wave-front, v = 0, hence 


U,aty ClU, Q 


wtere 

** 1=1 

From equation (10), we have 


U, a = «* a 


( 10 ) 


( 11 ) 
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which gives dependence of u a at the wave front on time Since at wave front, t = \j /, 
u *, a is value of u, a at t = t* 


For unconstrained equilibrium case, 




f al M 

V a *o yj 


The critical time for shock-formation is given by 


t _^,2 1*\ (r 0 + i)g * / * ; 

o l(r o +l)ff*/*+oJ 

Velocity gradient is given by 


u * e~ 

W-P) 

l + r ° +1 u*, ;c 


2 

< a ; 


( 12 ) 


which m dimensionless form can be written as 


s = {^ + i( e -.,|" (13, 

For rarefraction wave-front, velocity-gradient at an instant /* is positive and will 
decrease exponentially and damped out ultimately This nature is exhibited by curve 
d and e m Fig. 1 For compressive-wave front, velocity gradient at the wave front 
will be steepen up into a shock-wave as shown by curves ‘C and V in Fig 1 However, 
if 



CO 
VI 

Si 

0 

« 

> 
v> 
c 

s. 

x 
Ul 

1 

^ ^-Growth and decay of non-linear waves with planar symmetry 
u, x = 

r 0 +i 

velocity gradient at the wave-front will remain constant in time. 

( 11 ) Particular cases 

Case I - for « = 0, we have, law of propagation for ordinary gas 

Case II - for n =1, we have law of propagation for single relaxation mode 
discussed by Clark 9 and Rarity 10 

Case III - for n = 2, case for disociation and vibrational mode can be discussed. 

(in) Cylindrical and spherical wave-front 

For cylindrical and spherical wave front, we put v=l and 2 in equation (9) 
respectively and accordingly have equations 



8 = [(1 + 2t|) 1/2 e cn {l + MF(ri)}]' 1 
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and 


8 = [(1 + 2rj)e OT1 {l + Xff (r))}]" 1 

respectively Nature of cylmdncal and spherical wave-fronts are almost similar to that 
for plane-wave front 


Weak Non-Linear Wave 

Here we discuss certain properties of weak-non linear wave in reacting-gases for 
planer-wave front and obtain exact solution with the help of Laplace-transformation. 

Following Chu 15 governing equations in this case are 



Applying Laplace-transformation and initial boundary conditions to equations (14) to 
(18) and on inversion, we have 

f: * 


( 19 ) 
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Cl ■C+' oa 

u(a,v)=V + ~ J/© e «-w^ + 0 (£ 2 ) ( 20 ) 

2TCI / — joo 


€ 

rv _L x i f _L_ — 

V 

*r ft) 

» al 

r+2+ S^,+i 

f 1 _ 1 N 
< a ‘‘ a f ; 


(XT \y i 

2%i 

J \ 

Z-CLoo A* 

,+ 6 m '4 

i n 

t J 

1 

I 


= 0(e 2 ) (21) 

where £ is a positive number which satisfies the condition that all the singularities of 
these integrals are to the left of line \- £ m complex ^-plane. Inverting equations 
(20) and (21), we have 

a = a (x,t) and vp = ip (x, t). 


which give the trajectories of an outgoing wave and particle path in (x,t) plane. 

Thus 

a = / - x+ e a (!) (x, t) + 0(e 2 )| ^ 

\j/ = r+ei|/ (1) (x,/) + 0(e 2 ) j 

Equations (19) to (21) describe the trajectories of an outgoing wave and particle 
paths in (x, t) plane The formation of shock is characterized by t, a = 0. 

Central Expansion Fan 

We consider the expansive action of a receding piston with 



U(o 


for t < 1/p 
for t >l/p 


(23) 
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where g(l/P) = -1/2 p, g' (1/p) = -1 and P is dimensionless piston acceleration. Here 
we consider a particular case in which the piston velocity changes from 0 at t = 0, to 
-e as P -> co 

Introducing a new characteristic level z defined by z = Pa, for a < 1/p, we see that 
as a vanes from 0 to 1/p, z will change from 0 to 1. As we are interested m a solution 
to the problem in which an abrupt change in piston velocity generates expansion- 
wave, following Chu 15 , we have 


„(» ] 

v a) _ 

Z > 

x$ = u il) +ay+t«\ 

(24) 

p,?-<-~p,y= o 

(25) 

?, m = C,, po p 0 \ i = 1,2, ,» 

(26) 

,=1 x, 

(27) 

4’ =-4-^-^} / = 1,2,...* 

r T z 

(28) 


(29) 


i=i 


Solutions of equations (24) to (29) in limit when p — > oo, can be written as 

x (1) = c, (1) = 0, z = l,. ... ,n ►, 

^lo+IJU -mA 

2 n v ' I 


( 30 ) 
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where T and Q have the same meaning as before Thus complete solution m 
expansion fan valid upto first order in e is 

p = Po -^e-™ + 0(e 2 ) 

u = -ePe“ n ' f +0(€ 2 ) 

S = S o +0(e 2 ) 

c, =c, 0 +0(e 2 ),; =1,2,.. ,n ( 31 ) 

x = v f / + 0(e 2 ) 

t = V+ - e- m ) + 0(e 2 ) 

2 O J 


The characteristic for this case is given as 

r 0 + iep 


t = x + - 


2 n. 


(l-e- UY ) + 0(e 2 ) 


( 32 ) 


For xO«l, i.e for time small compared to the relaxation time, we have 


x = /<|I + I^llep}> +0(e 2 ) 


For xQ»l, characteristic in expansion fan becomes 

r o +leP 


( 33 ) 


/ = x + - 


2 a 


+ 0(e 2 ). 


( 34 ) 


The velocity field in the physical space for this case is given by 


u = 


f \ 


v r ° + 1 ; 


Qe-^x-Ol 


l~e 


-Ox 


+ 0(e 2 ). 


( 35 ) 


Thus, velocity gradient at leading wave front x = t, is given by 
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2 [ Qe-°* } , 2 , 

“'■“'rwlPY^r 

In limit when x, -» oo, Q 0, and we have 


u(x,t) = 


f ne ^ 


r„ + i 


1-e 


-at 


+ 0(e 2 ) 


(36) 


(36) 
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Abstract 

Mathematical model analysis of turbulent free jet envisaging the physical processes like 
Brownian diffusion, drag force and volume fraction m relation to dust loading in the 
atmosphere has been studied for better monitoring and prediction The analysis presented the 
modified form of Reynolds equations for turbulent boundary layer to account for the 
interactions due to physical processes that arise due to presence of SPM Crank-Nicholson 
implicit finite difference scheme has been employed for model validation The model 
prediction depicts the interdependence of SPM and mixing length The analysis is m good 
agreement with Schhchting and experimental results of Reichardt Consideration of volume 
fraction and Brownian diffusion of SPM helps m the migration of SPM through a longer 
distance 

(Keywords turbulent free jet/suspended particulate matter) 

Introduction 

Flow of fluids with SPM has been the focus of the study of engineers and 
meteorologists for its wide range of applications. Most of the practical fluid flow 
problems involve turbulence. There are a number of laminar flow problems for which 
the equations of motion have been solved exactly and many more for which the 
equations can be solved by certain approximations without greatly affecting the 
validity of the results. However there does not seem to have exact solution for 
turbulent flow. The approximate equations describing turbulent flow depend on so 
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many assumptions that it is difficult to say whether agreement with experiment is the 
result of reasonable simplifications or of fortuitous cancellations of the errors arising 
from the assumptions. In spite of the difficulties of obtaining a complete theoretical 
solution of the turbulent flow problem, useful quantitative relations have been 
obtained by a combination of theoretical reasoning and empiricism. 

Kliegel 3 and Soo 4 have excluded turbulence in their studies of dynamics of fluid 
with SPM. Rose 5 , Soo 6 and Wen et al 1 have studied similar problems with 
empiricism. Soo 8 has studied the fully developed turbulent pipe flow by extending the 
l/7th velocity distributions of Schilchting 1 to the gas solid suspension. Soo 9 and 
Ryhming 10 have considered two dimensional jet mixings of dusty fluid with a clear 
fluid. While Soo 9 did not consider the diffusion of SPM in his analysis, neither Soo 9 
nor Ryhming 10 considered the Stokes drag terms in the fluid momentum equations. 

The present work proposes to study the turbulent boundary layer equations 
governing the jet flow with SPM in two different approaches. In the first approach, 
mixing length is assumed to be unaffected by the presence of SPM. In the 2nd 
approach, when SPM is assumed to affect the mixing length, the Reynolds stress 
terms are replaced by Boussinesq approximations. The interactions of SPM with the 
carrier fluid, which is summarized as those forces like drag-forces due to slip, 
transverse forces due to slip shear, lift forces due to Brownian motion of the SPM, 
have been considered. In addition to these forces, volume fraction of SPM, which 
plays an important role, has been considered. Brownian diffusion of SPM is also 
considered due to small size of particles. Prandtl’s boundary layer approximation for 
simplifying Reynolds equations for a turbulent boundary layer flow and Prandtl’s 
mixing length theory and Boussinesq approximations have been adopted to close the 
system. 


Mathematical Formulation 

Typically, the jet is turbulent and is characterized by a uniform velocity profile at 
the nozzle exit. However with increasing distance from the exit, momentum exchange 
between the jet and the ambient causes the free boundary of the jet to broaden and the 
potential core, within which the uniform exit velocity is retained, to contract. Down- 
stream of the potential core the velocity profile is non-uniform over the entire jet 
cross-section and the maximum velocity (center) decreases with increasing distance 
from the nozzle exit. The region of the flow over which conditions are unaffected by 
the impingement surface is termed as the free jet. 
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A realistic theory of the motion of fluid with SPM must account for (i) the solid 
fluid mixture is heterogeneous, (si) the average velocity of the particles near the wall 
is not zero and the excess of particle velocity over the fluid velocity is positive near 
the wall, ( 111 ) the concentration of the solid particles is not uniform because the solid 
particles have greater density than fluid 

Basically there have been two approaches to the problem of fluid-solid flow, viz. 
(i) extending the dynamics of a single particle to the case in a similar manner as is 
done in case of molecular theory, (ii) modifying the continuum approach to account 
for the particles in the fluid phase. In the second approach, the dispersed solid phase is 
viewed as a continuum despite the discreteness of the particles Following the second 
approach, a number of authors viz. Soo 4 , Marble 11 ' 12 , Saffinan 13 , Batchelor 14 , Hinze 15 
and others, have presented modified forms of N. S. equations for the fluid and particle 
phases. 

In case of fluid velocity ranging from 15 to 30 meters per second, diameter of 
solid particles ranging from 100 to 200 microns (for glass beads) and solid to air mass 
ratios of 0.05 to 0.15, compressibility of fluid phase and the effect of gravity on the 
density distribution of solid particles are negligible and velocity distribution of fluid 
stream is not significantly affected by solid particles. 

The equation of continuity for the pseudo fluid of SPM can be written as 
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^ , p *“ J + f ( +P * V ') =0 

or £^) + |;(p'p v i>) =0 


since is constant, the equation of continuity reduces to 


du r 


- + - 


cv 




or 




( 2 ) 


The diffusion of the species of pseudo fluid of SPM will be governed by law of 
conservation of species and the Fick’s Law can be used to determine the species 
diffusion rate at any point in the medium and the resulting equation could be solved 
for species concentration distribution. The species conservation equation is given by 


dc 5c 

u — + v — 
dx dy 


= D r 




f &c_ &c 
dx 2 + A - 2 




( 3 ) 


where c = — — ,p -density of particle phase, p m -density of the mixture of fluid and 
P m 

particles 


so 


Pm 


z Pr+P\ 



Pd 

where — = <j) is the volume fraction. 

Pjrp 

Considering finite volume fraction and diffusion of SPM through carrier fluid, the 
governing Reynolds equations in two- dimensional cartesian coordinate system are 
given by 
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^dc +y <3c _ 3 

dx 8y dy 



( 8 ) 


where the variables u, v, u p , v p , p, p, p p , c etc represent the mean quantities. 

In deriving the above equations the following considerations are accounted for 

The terms (u p -w'jand (v^ -v ) are negligible in comparison with (u p -uj and 

Since, by the time the flow becomes turbulent the SPM have had enough time to 
cover many times of the particle velocity relaxation length so that the slip velocity 
(u p - u ) and (v p -v) become very small and the perturbation quantities u ,u p ,v ,v p 
are such that we can assume u p =o[u'\v p =o(v j 

As has been done in denvmg boundary layer equations, where 
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8_ 

dx 


0 



rn 

, 5 , 


The term, — ( v' 2 - (l - <J>)r' 2 ] , can also be neglected 
dx\ ) 


and 8) 

P dy 


1 dp ' -■ g 1 °£ 
p dy ay ^ ' p St 



If we assume that the SPM does not affect the mixing length as the particle 


concentration is very lean, terms 


i(-^ 


and 


8 _ 

dy 


can be neglected. 


Using Prandtl’s nuxing length theory the Reynolds stress, -u'v , is replaced by 


r 


^du } 2 

dy 


where / is known as the mixing length which is the distance up to which 


migration of fluid elements is assumed to take place. 


By analogy with Fick’s Law of diffusion, we write 


vc 



in which K c is the turbulent diffusivity or ‘eddy diffiisivity’. 

After making the above mentioned replacements for the turbulent momentum flux 
and turbulent mass flux and introducing non dimensional variables 
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The boundary conditions are 

«(0,y) = l u p (0,y)=l 

u(x, oo) = 0 w JJ (x,oo) = 0 


w y (x,0)=0 
w (x,0)=° 


the governing equations (9)-(13) are solved by using Crank-Nicholson finite 
difference scheme 19 . 
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However, if we assume that the SPM dees affect the mixing length, the 2nd terms 


of right hand side of eqn (6) and (7) i.e 



can not be neglected 


Boussmesq, in 1877, had suggested that the turbulent stress in analogy to viscous 
stress can be related to the rate of strain, in case of incompressible flow, as 


— u,u , 


:K m 


' du, t du^ 

to' , 


where K m is the eddy diffusivity for momentum transfer for fluid phase. Boussnesq 
approximations for particle phase may be stated as follows 


■ U P, U P ] = K mp 


du 


Pi 


Sr, 



to, 


J 


where K mp is the eddy diffusivity for momentum transfer for particle phase 

Introducing non dimensional variables, as earlier, the eqn. (4) to (8) can be written 
as(after dropping the stars) 


du dv . 
— + — = 0 
dx dy 


M) 

♦ 


du du 

u — + v — 
dx dy 


= (l + (l-$)s m Re)^ + ap,(», - u ) 


top , du 

U n ~ + V„ ~ 

p dx p dy 


= 4,-«)+(K* Re-sfe- 

dy 


4> 


top , dv 
u . — - + v„ — — 
p dx r dy 


= 4, - V )+( 2 K, Re- 

dy 


(15) 

(16) 

(17) 


(18) 
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d Pj> 


d p p „- d 2 p f 


u — — + v—-^- — s' c Re — -y- 
dx dy 8y 2 


V_ K m K 

where 8 = — £„ f — , s„, 7 . — • „ . 

v ” Uk mp Uk 


mp 


tP,+K.) 


Uk 


(19) 


and are solved by using the Crank-Nicholson scheme. 

Discussion of Results 

The numerical results obtained suggest some of the characteristic featuies of gas 
particulate flow in the turbulent jet The values used for numerical computations are 
presented m Table-2 which have been modified from the data used by Fenton and 
Stukel 16 , given in Table- 1, using the following relations, to suit the present problem. 



X - 2 

M 

a 2 U 

K m d 


K K m „ 


K c „ 


9 

Ip J 

V 

~ Ud V 

C 

Uk~ K 

m 


V ° m 

A. w 

tn 

Table 1- 

Data used by Fenton and Stukel 16 





a 

d 

(m) 


ft , 

(Kgirr) 

a 

(micron) 

u 

fins'*) 

K 

m_ 

Ud 

K 

mp 

Ud 

K 

Ud 

0 0141 

0 003175 

801 

5 Q0E-05 

12 68 

0 0481 

0 489 

0 632 

0 0141 

0 003175 

1602 

5 00E-05 

12 68 

0.0481 

0 489 

0 632 

0 00657 

0 003175 

801 

5 00E-05 

18.07 

0.0481 

0.469 

0 675 

0 00657 

0 003175 

1602 

5 00E-05 

18.07 

0.0481 

0 469 

0.675 

0 0157 

0 003175 

801 

5 00E-05 

18 26 

0.0443 

0 551 

0.686 

0 0157 

0 003175 

1602 

5 00E-05 

18 26 

0.0443 

0 551 

0 686 

0.0146 

0.00635 

801 

5.00E-05 

27 

0.0418 

0 608 

0.804 

0 0146 

0 00635 

1602 

5 00E-O5 

27 

0 0418 

0 608 

0 804 
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Table 2- Data used in the present analysis 


a 

p 

(Kgm' 3 ) 

p 

(Kgm J ) 

p 

(Kgm- 3 ) 

a 

(m) 

u 

(ms -1 ) 

X 

(m) 

K 

rn_ _ 

UK m 

K 

mp 

UK mp 

K 

£_ _ g 

UK " 

0 0141 

0 9752 

1 54E-05 

89: 

1 00E-04 

12 68 

1 46E+0Q 

1 04E-04 

5 1 0E-05 

6 59E-05 

00141 

0 9752 

1 54 E-05 

1602 

1 00E-04 

12 68 

2 93E+00 

5 22E-05 

2 55E-05 

3 30E-05 

0 00657 

0 9752 

1 54E-05 

801 

1 O0E-04 

18 07 

2 09E+00 

7 32E-05 

3 43E-05 

4 94 E-05 

0 00657 

0 9752 

1 54E-05 

8010 

1 Q0E-Q4 

18 07 

2 09E+01 

7 32E-06 

3 43E»06 

4 94E-06 

0 0157 

0 9752 

1 54E-Q5 

801 

1 00E-04 

18 26 

2 11 E+0G 

6 67E-05 

3 68E-05 

4 58E-05 

0 0157 

0 9752 

1 54E-Q5 

1602 

1 Q0E-Q4 

18 26 

4 22E+0Q 

3 34E-05 

1 84 E-05 

2 29E-05 

0 0146 

0 9752 

1 54E-05 

801 

1 OOE-04 

27 

3 12E+00 

8 51 E-05 

5 18E-05 

6 85E-05 

0 0146 

0 9752 

1 54E-05 

1602 

1 00E-04 

27 

6 24E+00 

4 26E-05 

2 59E-05 

3 42E-05 


Fig 2 compares velocity profile of carrier fluid, which is numerically computed 
and drawn against the velocity profile of Schlichting', viz. 


1 ~ tar) h 2 v - 2ri(l - tanh 2 r|)~ 



The computed profile agrees well in the core of the jet and differs from the profile 
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Fig 3 presents the velocity profile of SPM drawn against the mean velocity 
distribution of solid particles approximated by an equation analogous to the velocity 
profile of fluid The computed results show that the particle velocity decreases 
monotomcally from the value unity at the center of the jet to its asymptotic value in 
the ffee-stream 



Fig 4 shows the computed density profile of SPM drawn against that akin to the 
postulated velocity profile of fluid. The density profile agrees well in the core of the 
jet but deviates outside the core 



Fig 4- Density profile of SPM 
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Fig. 5 presents the velocity profile of carrier fluid in presence of SPM, which 
affects the mixing length. The Reynolds stress terms appearing in the governing 
eq uati ons have been approximated by quantities as predicted by Boussinesq. It is 
observed that when SPM affects the mixing length, the velocity profile of fluid 
deviates from the one, which have been computed using mixing length theory (Fig. 2), 
The fluid velocity falls sharply, from the maximum value of unity at the center of jet, 
inside the core of the jet and afterwards maintains a steady decline outside the core 



Fig. 5 - Velocity profile of fluid in presence of SPM when SPM affects mixing length 

Fig. 6 makes a comparison of the velocity distribution of SPM and that of the 
postulated distribution when the mixing length is affected by SPM. In this case also, 
the velocity distribution differs from the one that have been computed using mixing 
length theory. The velocity fells monotonically from the maximum value at the center 
of jet to zero in the free-stream. 



Fig. 6- Velocity profile of SPM when SPM affects mixing length 
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Fig. 7 depicts the density profile of SPM which when compared with that obtained 
in the 1st approach (Fig. 3) resembles well. No significant change is observed in the 
density profile in the two approaches. 



Conclusions 

A comparative study of the results obtained in both the approaches suggests that 
the first one, where mixing length theory has been used, gives a better structure of the 
boundary layer flow. 

The deviations marked in the computed results from the experimental results 
might have been caused by the following factors 

• Boundary layer assumptions 

• Postulation of velocity and density distribution from experimental results 

• Assumptions made to close the system 

The shortcomings of the present result may be overcome by considering higher 
order closure schemes. 
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Abstract 


The Rayleigh-Taylor instability of a Newtonian viscous fluid overlying Walters' (model B') 
viscoelastic fluid is considered For the stable configuration, the system is found to be stable or 
k 

unstable if v' < or > — — The system is found to be unstable in the potentially unstable 

i ea 

i 

case. The variable horizontal magnetic field is also considered For the stable configuration, in 

k 

the hydromagnetic case also, the system is found to be stable or unstable if v' < or > — — For 

i sa 

i 

the unstable configuration, the system is unstable for the wave-number band k < k* or if vj > 
k 

— — The magnetic field has stabilizing effect and completely stabilizes the wave number band 
i 

k 

k> k if v' < — — . This is in contrast to the stability of Newtonian superposed fluids in 
1 ea 

i 

porous medium, where the system is always stable for the stable configuration. 

(Keywords * viscoelastic (Walters' B' )-Newtonian fluids/variable magnetic field/ porous 
medium) 
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Introduction 


The stability derived from the character of the equilibrium of an incompressible 
heavy fluid of variable density (i.e. of a heterogeneous fluid) is termed as Rayleigh- 
Taylor instability. The Rayleigh-Taylor instability of a Newtonian fluid under varying 
assumptions of hydrodynamics and hydromagnetics has been studied by several 
authors, and Chandrasekhar 1 , in his celebrated monograph, has given a detailed 
account of these investigations. The Rayleigh- Taylor instability problems arise in 
oceanography, limnology and engineering. The problem of Rayleigh-Taylor 
instability in fluids in a porous medium is of great importance in geophysics, soil 
sciences, ground water hydrology and astrophysics. 


With the growing importance of non- Newtonian fluids in modern technology and 
industries, the investigations on such fluids are desirable. There are many elastico- 
viscous fluids that cannot be characterized by Maxwell's constitutive relations or 
Oldroyd's constitutive relations. One such class of elastico-viscous fluids is Walters' 
(model B' ) fluid. In recent years, the investigation of flow of fluids through porous 
media has become an important topic due to the recovery of crude oil from the pores 
of reservoir rocks. The flow through porous media is of considerable interest for 
petroleum engineers and for geophysical fluid dynamicists. A great number of 
applications in geophysics may be found in a book by Phillips 2 . When the fluid 
permeates a porous material, the gross effect is represented by the Darcy's law. As a 
result of this macroscopic law, the usual viscous and viscoelastic terms in the equation 
of Walters' (model B') fluid motion are replaced by the resistance term 



where p and p' are the viscosity and viscoelasticity of the 


Walters' (model B') fluid, k\ is the medium permeability and u is the Darcian (filter) 
velocity of the fluid. Sharma et a /. 3,4 have studied the thermosolutal convection in 
Walters' (model B') fluid in porous medium in the presence of magnetic field and 
rotation, respectively. 


Generally, it is accepted that comets consist of a dusty 'snowball' of a mixture of 
frozen gases which, in the process of their journey, changes from solid to gas and 
vice-versa. The physical properties of comets, meteorites and interplanetary dust 
strongly suggest the importance of porosity in astrophysical context (McDonnell 5 ). 
Wooding 6 has considered the Rayleigh instability of a thermal boundary layer in flow 
through a porous medium. In stellar interiors and atmospheres, the magnetic field may 
be (and quite often is) variable (and non-uniform) and may altogether alter the nature 
of the instability. Sharma and Sunil 7 have studied the Rayleigh-Taylor instability of a 
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partially ionized plasma in a porous medium in presence of a variable magnetic field. 
Sunil and Chand 8 have studied the Rayleigh-Taylor instability of a plasma in the 
presence of a variable magnetic field and suspended particles in porous medium. 
There is growing importance of non-Newtonian fluids in chemical technology, 
industry and geophysical fluid dynamics. 

Keeping in mind the importance of non-Newtonian fluids in modem technology 
and various applications mentioned above, the present paper is devoted to the 
consideration of the stability of superposed Newtonian-Walters' (model B') fluids in 
porous medium. The effect of variable horizontal magnetic field, bearing relevancy in 
geophysics, is also considered separately. 

Formulation of the Problem and Perturbation Equations 

Consider a static state in which an incompressible Walters' (model B') fluid layer 
of variable density is arranged in horizontal strata and the pressure p and density p are 
functions of the vertical coordinate z only. The character of the equilibrium of this 
initial state is determined by supposing that the system is slightly disturbed and then 
by following its further evolution. The fluid is under the action of gravity g (0, 0, - g). 
This fluid layer is assumed to be flowing through an isotropic and homogeneous 
porous medium of porosity e and medium permeability k\, where the porosity is 
defined as 

e = (volume of the voids)/ (total volume), (0 < s < 1). 

For very fluffy foam materials, e is nearly one and in beds of packed spheres, e is in 
the range 0.25-0.50. 

Let p, p, pi, p' and it (u,v,w) denote, respectively, the pressure, density, viscosity, 
viscoelasticity and filter velocity of the pure Walters' (model B') fluid. Here g is the 
acceleration due to gravity, x- (x, y, z ) and X = (0,0,1). Then the equations of 
motion and continuity for the Walters' (model B') fluid are 


P 

£ 




v 1 


-Vp-pgX- 


p-p 


, 8 _ 

dt 


u, 


( 1 ) 


V.5 = 0 


( 2 ) 


Since the density of a fluid particle remains unchanged as we follow it with its 
motion, we have 
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s — + (u v)p = 0 
dt v ^ 

( 3 ) 

Consider a small perturbation on the steady state solution, and let Sp, Sp and u (u, 
v, w) denote, respectively, the perturbations in pressure p, density p and fluid velocity 
(0,0,0) Then the lineanzed perturbation equations governing the motion of the 
Walters' (model B') fluid layer through porous medium are 

- ~ = -vfy ~ g 5 P% ~ “f V ~ & J V, 

s dt k x \ dt ) 

( 4 ) 

< 

1! 

o 

( 5 ) 

ps 

s— 8p = -w(Dp) 

(6) 

Dispersion Relation 


Analyzing the perturbations into normal modes, we assume 
quantities have an x, y and t dependence of the form 

that the perturbation 

exp(i^x + ik y y + nt\ 

( 7 ) 


where k x , k y are the wave numbers along the x- and y- directions respectively, k = 
■J(k x + k y ) is the resultant wave number of disturbance and n is the growth rate which 
is, in general, a complex constant. 

For perturbations of the form (7), eqns (4)-(6) give 


n 1 . . . 

— + — - (v - v n) pu = 5o, 

s k x 
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n 1 , , . 

— i (v-vn) 

s 


pw = -D&p - gSp, 


( 10 ) 


tk x u + iky v + Dw - 0, (11) 

ew5p= -w(£>p), (12) 

where v = — , v' = — and D = — . 
p p dz 

Multiplying (8) by -ik x , (9) by -ik y , adding and using (1 1), we obtain 


n 1 . , . 

— + — (v-vn) 
s /c. 


p Dw = -k 2 hp. 


(13) 


Eliminating 8 p between (10) and (13) and using (12), we get 
—[D(pDw) - t 2 pw]+ -i- - \i!n)Dw} - k 2 (p - p/«)w]+ (Dp)w = 0 (14) 

8 /C| S Yl 

Two Uniform Newtonian and Viscoelastic (Walters' B') Fluids 
Separated by a Horizontal Boundary 

We consider the case of two uniform fluids of densities, viscosities, 
viscoelasticities; p 2 , p. 2 , 0 (upper newtonian fluid) and pi, p,, pj (lower Walters' B' 
viscoelastic fluid) separated by a horizontal boundary at z=0. Then, in each region of 
constant p, constant p and constant p', eqn. (14) becomes 

(D 2 - k 2 )w =0 (15) 

The general solution of eqn. (15) is 

w = Ae +k2 + Be kz , 

where A and B are arbitrary constants. 


( 16 ) 
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The boundary conditions to be satisfied in the present problem are : 

( I ) The velocity w should vanish when z -> + oo (for the upper fluid) and z - 
oo (for the lower fluid) 

( II ) w(z) is continuous at z = 0 

(ui) The jumps condition at the interface z = 0 between the fluids is obtained by 
integrating eqn (14) across the interface at z =0 and is 

n \ k 2 

~[P 2 Dw 2 -Pl^lLo “(^1 -K«)^lL=0 = “ — (P2 -PlK» (17) 

E ATj E/1 

remembermg that upper fluid is Newtonian and lower fluid is Walters B' viscoelastic 
Here w 0 is the common value of w at z =0 

Applying the boundary conditions (l) and (ii), we can write 


W, =Ae +k \ (z < 0), 

(18) 

w 2 = Ae~ k \ (z.> 0), 

(19) 


where the same constants has been chosen to ensure the continuity of w at z =0 
Applying the condition (17) to the solutions (18) and (19), we obtain 



K 


I n 2 +- S 


k \ 


(v 2 a 2 + v t a, )n - gk( a 2 - a, ) = 0 , 


where 


tt l,2 _ 


Pi, 2 

Pi +P2 ’ 




( 20 ) 


(a) Stable case (p 2 < pi) 

For the potentially stable arrangement (p 2 < pO, the system is stable or unstable 
according as 
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v{ < or > 


A. 


( 21 ) 


(b) Unstable case (p 2 > Pi) 


For the potentially unstable arrangement (p 2 > pi), there is atleast one change of 
sign in eqn (20) and so this equation has one positive root. The occurrence of positive 
root implies that the system is unstable. 

Effect of Variable Horizontal Magnetic Field 

Here the problem and configuration is same as described above except that the 
incompressible Walters' (model B') fluid layer arranged in horizontal strata is acted on 
by a variable horizontal magnetic field H (Ho(z), 0.0). Thep the equations of motion 
and the Maxwell’s equations are 


P 

e 




-Vp~pgl-~\ p-p' 


K 


\\ 


dt 


a + — (v x )x H, 

4tc ' ’ 


( 22 ) 


V.# = 0, (23) 

s ^ = {S.v)u-(u.V)H, (24) 


together with eqns. (2) and (3). here \x e denotes the magnetic permeability. 

Let h (h x , h y ^ h z ) denotes the perturbation in the variable horizontal magnetic field 
H (H q (z), 0.0). Then the linearized perturbation equations, governing the motion of 
the Walters' (model B') fluid layer through porous medium in the presence of variable 
horizontal magnetic field are 


p du „„ i 

7 ¥ = -V6p-g8 P A,-- 


, a 

v dtj 


!m + -^[(Vxa)x^ + (Vx^)xa], (25) 

4n 


vi = o, 


(26) 
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s- = (Bv)a-(n.v)R, 


together with eqns (5) and (6) 

For perturbations of the form (7) eqns. (25)-(27) give 


7 + f (v ~ V '”) P w = ~ ,k x&P + ~~h z (DH 0 ), 


■ + -(v~ v'n) pv = - lky 5p + Ml 0 k x h y - Jky h x ), 


• + — (v - v'n) pw = -DSp + M±\, ,k x h z - Dh x -h j^ 


h „ ——1 _ 


ik x h x + / kyhy + Dh, = 0, 
e«Ax = ik x H 0 u-wDH 0 , 


znhy — ik x Hov, 


g5p, ( 30 ) 


together with eqns. (1 1) and (12). 

Equation (29), with the help of (32) and (33), becomes 


- + -(v-v'n)J \. pv = - tky 8p + ^ (lW + J k yW DH 0 ), 


where C, ik x v- i ky u is the z-component of vorticity. 
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Multiplying (28) by -ik x , (35) by -ik y , adding and using (1 1), we obtain 


" + l(v- Vn)]pDw = -k 2 hp + (DHo )w _ h , (dH 0 ) 

s k\ 4tis« 4ns n 4 n 


Eliminating Sp between (30) and (36) and using (1 1), (12) and (28)-(34), we get 


D(pDw)~ k 2 pw ]+ 7“ M(P “ l 1 '”)} - k 2 (p - p'«)> 


d(h 2 Dw)~ H 2 k 2 w\-^~ (Dp)w 


Here also we consider the case of two uniform Newtonian and viscoelastic 
Walters' (model B') fluids separated by a horizontal boundary in hydromagnetics 
exactly as in previous section. 

The jumps condition at the interface z = 0 between the fluids is obtained by 
integrating eqn. (37) across the interface atz = 0 and is 

-\p 2 Dw 2 -p l Dw x \ =0 +~[x 2 Dw 2 -(pj -P^IMLo 

£ /Ci 


— (P 2 -PlK 

4tts n s n 


Following the procedure as in previous section, we obtain 


1 n 2 +y-(v 2 a 2 +v l a. l )n + \k 2 x Vl - gk{a 2 -aj^O 


where for the sake of simplicity, we have considered that the Alfven velocities of the 
two fluids are the same, so that 
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V 2 = 

v A 


yjjj = Mg 

47ip, 47tp, 


(a) Stable case (p 2 < pi) 

For the potentially stable arrangement (p 2 < pi), the system is stable or unstable 
according as 


Vj < or > 


ECC, 


(40) 


Thus, for two uniform Newtonian and viscoelastic Walters' (model B') fluids in 
porous medium in hydromagnetics and for the potentially stable case, the system is 

k 

stable or unstable according as v( < or > — — 

ea, 

This is m contrast to the stability of Newtonian superposed fluids in porous 
medium, where the system is always stable for the stable configuration 


(b) Unstable case (p 2 >pi) 


For the potentially unstable arrangement (p 2 > pi), if 
k 2 V] > gk{ a 2 - a,) and v{ < 

sa, 


(41) 


Eqn (39) does not admit any change of sign and so has no positive root. Therefore, 
the system is stable However, if 


kxVj <gk(a 2 -a x ), (42) 

the constant term m eqn (39) is negative. Eqn. (39), therefore, allows one change in 
sign and so has one positive root. The occurrence of positive root implies that the 
system is unstable 
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k 

Thus for the unstable case (p 2 > pi) and v[ < — — , the magnetic field has got 

SC6] 

stabilizing effect and completely stabilize the system for all wave numbers which 
satisfy the inequality 


k\v] >gk(a 2 - a,). 


(43) 


i.e. k> k* 


where k* = 


g(a 2 ~ a ;) 

y\ 


sec 2 0, 


and 0 is the angle between k and H . 


However, for the potentially unstable arrangement (p 2 > pi), the system is 
unstable for the wave-number band k < k* or if v[ > . 

SO-! 


Thus, for the potentially unstable configuration and for a uniform Newtonian 
fluid overlying a viscoelastic Walters' (model B') fluid in porous medium in 
hydrodynamics, the system is unstable for all wave numbers, whereas in 
hydromagnetics, the magnetic field has stabilizing effect and completely stabilizes the 


wave number band k > k*, where k* = &( a2 


V} 


sec 2 0 
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Abstract 

We study the problem of optimal implementation of parallel divide-and-conquer algorithm on 
binary de Bruijn networks. A divide and conquer algorithm is modeled as a temporal complete 
binary tree computation structure An important contraction property between two successive 
binary de Bruijn network is revealed A twice-size compute binary is mapped to a de Bruijn 
network. Two nodes in the compute binary tree are mapped to a single node The rnappmg is of 
dilation one, communication contention free and of good load balance 

(Keywords parallel divide/ conquer algorithm de Bruijn network/ parallel algorithms) 

Introduction 

Divide and Conquer method is an important parallel programming paradigm 1 ’ 3 . 
The problem of implementation of such algorithms onto various parallel architectures 
such as hypercubes, tree machines, mesh connected computer, etc. has been studied. A 
de Bruijn network 1 consists of n=2 k nodes. Let a k . i, a k . 2 , . . a h a 0 be the address of a 
node in the de Bruijn network The two nodes reachable via directed edges from the 
node are 


aw, • • •> ct \ , a 0 0 

a k~2, #0 1. 
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The diameter of a de Bruijn network 1 with 2 k nodes is k, which is about half the 
diameter of a shuffle-exchange network with the same number of nodes As an 
interconnection network, a de Bruijn network has many superior features such as 
bounded degree, logarithmic diameter and high connectivity Recently, there has been 
increasing interest in exploring such networks for multiprocessor system 2 . 

In this paper, we model a divide and conquer algorithm as a complete binary tree. 
To improve the efficiency, we map two non conflicting nodes in a complete binary 
tree to a single node in the binary de Bruijn network. We present a mapping which 
achieves dilation one and communication contention free. This is accomplished by 
utilizing an important contraction property of binary de Bruijn networks. 


The Model of Divide and Conquer 

A classic task graph structure for a degree-2 divide and conquer algorithm is a 
complete binary tree 4 where, the root first starts to divide the work into two 
approximately equal parts to its two children and once an interior node receives the 
work, it also divides it evenly to its children. After the leaves receives their work, they 
start to do compute and the results are aggregated level by level back to the root. We 
observe a very regular temporal activation pattern in such a computation : nodes and 
communication edges are active level by level . 

Formally, a complete binary tree with 2”' 1 nodes is denoted as CBT(«). We call 
the communication in the dividing frcm level i to z+1 (z=l,...,n-l) i-th communication 
phase and cell edges from nodes m level z to nodes in level z+1 edges in level z to 
nodes in level z+1 edges in level i. 

Four metrics are used to evaluate a mapping For the computation. Load 
Balancing is a measure for the processors load, defined as the maximum number of 
tasks which are mapped to a single processor. Processor Contention measures the 
degree of the contention for a processor among tasks We call a mapping Processor 
Contention Free if all tasks in a processor are active at different times. For the 
communication, Maximum Dilation 5 is defined as the maximum length of the paths of 
the de Bruijn net chosen for the communication edges in the CBT. Communication 
Contention is used to measure the contention among message traffic. We call a 
mapping Communication Contention Free if paths of the de Bruijn network chosen 
for edges active at the same time are unoverlapped 4 . 
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The Binary de Bruijn Network 

Definition . A binary de Bruijn network (graph) BDG(£) 6 is an undirected graph 
with 2 k nodes Each node is represented as a A-digit binary number a kA a k . 2 . . a 0 where 
a, e {0, r- 1}. Every node a kA a k . 2 a 0 is connected to the nodes : a k . 2 . . aye (by 
left shifting ) and xa k .\ a k . 2 ... a x (by right shifting ) where x s {0, 1 } . 

Property 1: A CBT(£) is a subgraph of a BDG(&). 

Proof We construct a CBT as follows start with 10*' 1 as the root of the CBT(/t), 
right shift to 010*' 2 to construct its left child, and right shift to 1 1 Q k ~ 2 to construct its 
right child Inductively , any node zlO m , where m> 0 and z is an arbitrary string of 
length k-m- 1, has left child 0zl0 m l and right child lzlO m '\ This procedure can be 
repeated until we reach a string whose rightmost bit is 1, which is the leaf of the CBT 
to be constructed. 

Nodes m the CBT (n+1) are labeled based on their labels in BDG (n+1). We call 
such a labeling of a complete binary tree de Bruijn labeling 7 . 

We now reveal an important contraction existing between BDG («) and BDG 
(«-l) 

We define a function F as follows : 

F(x i x n . x x n ) = x\ (1) 

where x, s {0,1} and x' = (x, © x,+i) for i < n. (2) 


Lemma L F is a two -to -one ... notion which maps both x x .. x„.\X„ and its binary 
complement strings x, x„ to the same string . 

Proof ■ It is easy to see that F maps both Xi . . . x„. x x n and x, .. .x n _, x n to the same 
string 

x'j x’ 2 . . ,xf { where x[ = (xr, © x,+i) for i < n because x ©y = x © y ■ Conversely, 
the strings X] x„_ix„ which are mapped to y x ... >vi are the solutions to the system of 
equations 
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Abstract 

1729 is associated with the famous Indian Mathematician Ramanujan This interesting number 
happened to be the number of the taxicab in which Professor Hardy visited him in the hospital 
and was astonished at Ramanujan's remark that 1729 was the smallest natural number 
expressible as sum of two cubes in two distinct ways involving relatively prime pairs 

(1 ) A 3 + B y = C 3 + D 3 = N with {AJB)= 1 ,(C,D)=1, I 3 + 1 2 3 = 9 3 + 1 0 3 = 1 729 

We propose to investigate (2) x 3 + y* + r 1 + « 3 = 0 for nonzero reals Through this we run into 
(1) No mathematical proof that 1729 is indeed the smallest number of the form (1) is claimed 

Introduction 

We introduce the following notations : 

A 

<x, y> = x 3 + y 3 a = a\ 02 ab = a x b 2 - a 2 b 1 


(«, m) = 1 the natural numbers n, m have no common divisors. 
] m, w [ open interval on the real axis. 


5, = - 0.70326 


R =] - 00 , + 00 [ 


/=]-!, Sd 
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S 2 =+2.36992 
Si = -1 42195 
S 4 = + 0.42195 
A* = ]-co, [ 

£*=}£, + «>[ 


C* =] - 00, Si [ 
D* =] S 4 , + » [ 
E* =] S U S A [ 

P * =]S 3 , 5,[ 
Q*=]S*,S 2 [ 


J=]S 4 ,1 [ 


Section 1 

If three real non-zero variables retain the same sign and one the opposite sign (2) 
takes the form 

(3) A 3 + 5 3 + C 3 = D 3 

Otherwise when two are of the same sign (2) becomes 

(4) A 3 + £ 3 = C 3 + D 3 

with A, B,C,D positive real numbers 

Section 2 


Let (x p y Jt Z } , U } ) j - 1,2 be arbitrary solution sets of (2) 

(5) Then (x, y, z, u) with x = xj + Xx 2 , etc can be a solution set of (2) provided. 

(6 ) x _ - xfo + y\y 2 + zfo + u?u 7 

K 

K* 0 with K = x x x\ + y x y\ + z x z\ + u x u\ 

Substituting X m (5) we have- x = EIK, y = F/K, z = G/K, u = HIK 

AAA 

(7 1) E=yxy + zxz + u xu 

V V V 
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AAA 

(7.2) F=zyz + uyu+xyx 

V V v 

A A A 

(7 3) G=uzu+xzx + yzy 

V V V 

A A A 

(7 4) H = x ux + y uy + z uz 

V w V 


(7 5) K = x x 2 + y y 2 + z z 2 + u u-, 

V v V v 

Section 3 

If K* 0, (E, F, G, H) will be a solution set for (2). (since <3,4> + < 5, -6> = 0, 
and <p, q> + <-p, -q> = 0) 

Choose (x h y u z u u{) = (3, 4, 5, -6) 

fe, yi, Z 2 , u 2 ) ~{p,q, -p rq) with p*0,q±0 

(8.1) E = 40 p 2 + 20 pq - 6q 2 

(8.2) F=32 p 2 + I6pq + 12 q 2 

(8.3) G = 24 p 2 - 20 pq - I0q 2 

(8.4) F[ — -48 p 2 -16 pq - %q 2 

(8.5) K=%p 2 -2q 2 

Let .S' = 2 /i/^ then: 

(8.1*) e = 5S 2 + 5S-3 
(8.2*) f=4S 2 + 4S+6 
(8.3*) * = 3S“-5S-5 
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(8 4*) h = -65 2 - 45 - 4 
(8 5*) k = 5 2 - 1 


{(x, y, z, u) = ( g , f, e, h) satisfies (2) if 5 * ± I } 
From g . the zeros of x are S\ = (5 - VS5 )/6 = - 0 70326 


and5> = 


£ ±^ - = + 2 36992 
6 


From e the zeros of z are S 3 = (-5 - V85")/10 = -1 .42195 and 

s 4 = + 0 42195. 

10 

Section 4 

(25 + l) 2 + 5 =/> 0 thus y is always positive. 

-6 (5 +1/3 ) 2 - 10/3 = h < 0 thus u is always negative. 

On C* U5*, also on E* three variables retain the same sign whereas one variable 
is of the opposite sign thus (2) takes the form (3) On p* U Q* two variables are of the 
same sign and the other two are of the opposite sign thus (2) takes the form (4). 

(-e/S\ -f/S\-g/S\ -hi S\) = (g*, h*, e*,f* ) 

Thus swapping S with / = 1/5 we still get solution set with the same constraints on 
t as on 5 The range - go < + oo can be contracted to -1 < S< + 1, since 5 ^ ± 1 Thus we 
can choose 5 m I U J in order to get the solution (2) in the form (4). By choosing 
rational S = p/q with (p, q) = 1 in I U J we can rid the equation (4) of fractions and get 
the solution as stipulated in (1) : A 3 + B 3 = C 3 + D 3 = N where A, B, C, D are distinct 
natural numbers with no common divisors. Fifty solutions are appended in Table 1 
Table 2 is actual verification 1 < N < 1729. 
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Table 1 


StFl 

S or t 

<AJB>= 

<C, D> = 


N 

A 

B 

- C 

D 

N 


1/2, -3/4 

1 

12 

9 

10 

1,729 

(1) 

2/3,-16/19 

23 

94 

63 

84 

842,751 

(2) 

4/5, -10/11 

35 

98 

59 

92 

984,067 

(3) 

10/17, -4/5 

23 

134 

95 

116 

2,418,271 

(4) 

5/7,-13/15 

51 

178 

115 

162 

5,772,403 

(5) 

3/4 

SI 

180 

113 

166 

6,017,193 

(6) 

13/20, -5/6 

45 

196 

133 

174 

7,620,661 

(7) 

6/13,-8/11 

9 

210 

161 

172 

9,261,729 

(8) 

3/5 

45 

246 

173 

214 

14,978,061 

(9) 

-5/7 

5 

254 

197 

206 

16,387,189 

(10) 

-6/7 

73 

270 

177 

244 

20,072,017 

(11) 

7/8,-17/18 

111 

268 

151 

258 

20,616,463 

(12) 

8/15,-10/13 

35 

298 

219 

252 

26,506,467 

(13) 

9/16, -11/14 

51 

348 

251 

298 

42,276,843 

(14) 

-7/8 

107 

356 

227 

326 

46,343,059 

(15) 

7/11 

89 

410 

282 

362 

69,625,969 

(16) 

3/7 

3 

4i4 

323 

334 

70,957,971 

(17) 

-7/9 

57 

430 

313 

366 

79,692,193 

(18) 

5/6 

167 

436 

255 

414 

87,539,319 

(19) 

-8/9 

147 

454 

283 

420 

96,753,187 

(20) 

-13/16 

89 

460 

321 

402 

98,040,969 

(21) 


Table 1 Contd.. 
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St=l 

S or t 


<A y B>= 

<C, D> = 


N 


A 

B 

C 

D 

N 


4/7 

73 

470 

337 

404 

104,212,017 

(22) 

12/19,-14/17 

111 

522 

359 

460 

143,604,279 

(23) 

10/11 

229 

522 

285 

508 

154,245,637 

(24) 

8/13 

111 

562 

391 

492 

178,871,959 

(25) 

-9/10 

193 

564 

345 

526 

186,595,201 

(26) 

6/7 

243 

606 

347 

580 

236,893,923 

(27) 

5/8 

133 

644 

445 

566 

269,442,621 

(28) 

-9/11 

133 

654 

453 

574 

282,078,901 

(29) 

4/9 

17 

694 

537 

564 

334,260,297 

(30) 

11/13 

271 

690 

399 

658 

348,411,511 

(31) 

-17/20 

189 

732 

485 

658 

398,974,437 

(32) 

5/9 

107 

706 

555 

654 

450,680,139 

(33) 

-11/12 

303 

820 

487 

774 

579,186,127 

(34) 

13/14 

389 

860 

461 

842 

694,919,869 

(35) 

11/16 

239 

908 

599 

818 

762,265,231 

(36) 

-11/13 

233 

926 

617 

830 

806,672,113 

(37) 

7/9 

317 

934 

573 

870 

846,635,517 

(38) 

-12/13 

367 

966 

567 

916 

950,859,559 

(39) 

5/11 

37 

1,046 

805 

854 

1,144,495,989 

(40) 

8/9 

437 

1,030 

573 

996 

1,176,180,453 

(41) 

7110 

295 

1,076 

703 

974 

1,271,439,351 

(42) 


Table 1 Contd 
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St=l 

S or t 

<A3>= 

<C, D> = 

N 


A 

B 

c 

D 

N 


13/17 

361 

1,098 

681 

1,018 

1,370,799,073 

(43) 

6/11 

147 

1,134 

827 

964 

1,461,450,627 

(44) 

-13/14 

437 

1,124 

653 

1,070 

1503,488,077 

(45) 

11/19 

189 

1,162 

829 

1,002 

1,575,734,797 

1461 

-11/15 

63 

1,174 

895 

966 

1,618,346,071 

(47) 

9/10 

555 

1,284 

707 

1,246 

2,287,828,179 

(48) 

-14/15 

513 

1,294 

745 

1,236 

2,301,725,881 

(49) 

16/17 

591 

1,282 

679 

1,260 

2,313,422,839 

(50) 


Table 2- List of natural numbers \< N < 1729 expressible as the sum of cubes of two distinct natural 
numbers 

For* = 2, 3,4, 5,6,7, 8, 9,10,11,12 

<1, *> = 9, 28, 65, 126,217, 344,513,730, 1,001, 1,332, 1 ,729, respectively 
For* = 3, 4, 5, 6, 7, 8, 9,10, 11 

< 2, *> = 35, 72, 133, 224, 351, 520, 737, 1008, 1339, respectively 


For* = 4, 5, 6, 7, 8, 9, 10, 11 

< 3, jc> = 91, 152, 243, 370, 539, 756, 1027, 1358, respectively 


For* =5, 6, 7, 8, 9,10,11 

< 4, * > = 189, 280, 407, 576, 793, 1064, 1395, respectively 
For * = 6, 7, 8, 9, 10, 11 

< 5, *> = 341, 468, 637, 854, 1 125, 1456, respectively 
For * = 7, 8, 9, 10, 11 

< 6, * > = 559, 728, 945, 1216, 1 547, respectively 


Table 2 Contd 
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Table 2 Contd 
For* = 8, 9, 10, 11 

<7,x> = 855, 1072, 1343, 1674, respectively 


For x = 9, 10 

< 8, x > = 1241, 1512, respectively 
Finally <9, 10 >= 1729 
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Abstract 

Type-2 or inverted multiple Dirichlet integral is extended for the real scalar variable case Then 
the matrix-variate analogues of the extended Dirichlet integrals, the associated extended 
Dirichlet models and extended Dirichlet densities are considered Various types of results are 
obtained which are mathematically and statistically interesting 

(Keywords : generalized Dirichlet model/real matrix-variate distributions/Jacobians of matrix 
transformation/type-2 model) 


Introduction 

The type-1 Dirichlet integral is given by 

i " L,<i, ocn + +J *<i x “'" 1 0 - ~ x k ) ak+rl dx\ a . . a dx k 


- r (cx ,)■ £(aj 

r(a,+...+a w ) 1 

fox 31 (a ; ) > 0, j = 1, £+1 where 9t(-) denotes the real part of (•). Hence the 
associated type-1 Dirichlet density is given by 


fM> x k )= ^+f+ a ^ x r\ x?-\\- Xx 


F(oci) ...(a k+] ) 


■Xk) 


ocAr+1 


-1 


( 2 ) 
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0 < x ; < 1,7 = 1, , k, 0 <x, + ,.+ xt< 1, 5R(o,)> O,; = 1,. , k+\, andT'ife, , x k ) =0, 
elsewhere The inverted Dmchlet or type-2 Dirichlet integral is given by 

r r ^- 1 +...+x,)- (ai+ +afc+i) = r(a,) r(a ^) (3) 

J ° J ° 1 k 1 k r(a,+ 3 -a^,) v ’ 


for 33 (a,)> 0,7 = 1 ,..., k + 1 , and the associated inverted Dmchlet or type -2 Dirichlet 
density is given by 




x k )= 


£Xgjjt +a i-+i) 

r(a,) r(a M ) 


,“1-1 r a k~ 1 


x k k d + x l + 


+ X,)' (a,+ + “* +1 > 


(4) 


for 0 < Xj < oo, 7 = 1 ,..., k and fi(x Iv . ..,x k ) = 0, elsewhere The matnx-variate 
extensions of type-1 and type-2 Dmchlet densities or the matrix-variate analogues of 
(2) and (4) are available in the literature, see for example Mathai 1 . They are the 
following : 


ri-y v’ "*"••• 4"0ti:+l) | v |ai 

-MX,- .,X k )~ — — — “A, 

r p (ai) ,r p (a k+i ) 


P+1 


‘••I* 


I “it 


.El i. 

^ /-x,-.. 


-X 


IOJt+1 — 


P+1 


*1 


(5) 


for 0 < X, = X' ] <I,0<Xi + . . +X k <1, 93 (Oj) > ----- ,7 = 1, ...,& + 1 and zero elsewhere, 
where the following notations are used. 

All the matrices appearing m this paper are real symmetric p x p matrices unless 
stated otherwise A prime denotes a transpose, |(-)| denotes the determinant of (•), tr(-) 
denotes the trace of (•), Xj = X' > 0 means XJ is real symmetric positive definite. 0< 
Xj< I means XJ = X' > 0 and I - XJ > 0 or that all the eigenvalues of XJ are in the open 
interval ( 0 , 1 ), and I denotes the identity matrix. The standard notation (a) means 
matrix-variate real gamma or 
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f l \ f n _ 1 A 

r,(a) = 7t 4 r(a)r a - - r(a - 1) .r a - , 91(a) > ^ 


'-1 


( 6 ) 


=L,oW a ^^ ,r( ^ 


( 7 ) 


where dX = dx t] wedge product of the — - differentials dx tJ 's. Matrix- 

variate version of type-2 density is given by the following • 

MXu- X k )= 


r p ( a i +- +a w) 

r p (a,). r p (a k+l ) 


w ai ~ 


p + 1 


p+\ 


..jx t p--|/+x 1+ ...+jir t f (ai++at+1) 


( 8 ) 


for Xj= X'j >0,y = 1, k, 91(0,) > - , j = 1 ,. .,/; + ! and zero elsewhere. 

Matrix-variate functions appear in a wide variety of problems in Statistics, 
Economics, Physics and so on When dealing with non-null distributions of likelihood 
ratio criteria m testmg statistical hypotheses on the parameters of one or more 
multivariate nonnal or Gaussian distributions, functions of matrix argument come in 
naturally. An example may be seen from Conradie and Troskie 2 . Applications of 
matrix-variate functions in various areas may be seen from Mathai 1 and die references 
therein Liouville models are generalizations of the scalar Dirichlet models. Matrix- 
variate analogues of Liouville models may be seen from Gupta and Richards 3 . 

The model that we are going to consider in this paper is a generalization of (8), 
along with its various properties. 

A Generalized Matrix-variate Type-2 Dirichlet Model 

Let Xj = Xj >0 bo p x p real symmetric positive definite matrices having a joint 
density which is a generalized real type-2 Dirichlet density of the following form : 
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f{Xu. ,X k )= 

\X k r^\l + X 2+ + X k \»'\l+X i+ .+X k \ h 

x.|/+X 1 + +Xk \ <ai++ak+ ^ +M - 


for % = X] >0, y=l, *, 91(0,) > 5R(o, +1 +. +o k+1 + p, + .+p t ) > £-1, 

/* z. 

7 = 1 , ,k and / (X 1: A).) = 0 elsewhere, where c k is the normalizing constant This 

normalizing constant can be evaluated by integrating out the matrices one by the 
starting with X\ Collecting the factors containing X, we have following- Let the 
integration over X\ be denoted by Lj. Then 

L x = £,_^ >o |X 1 p”~T|/ + Xi+ . + A r k p ( “ 1+ rat+1+pl+ + m dX x 


P + 1 


= |7 + X 2 + ...+ X,f (a,+ +a * +1+Pl+ J |X, | ai ~ 


!/ + (/+ X 2 + ... ■ + X k )"* JT, (/ + X 2 + . + X k )‘2 


-(“] + +a k+ \ +P] + +p fc ) 


clX, 


where (I-f-X \ - . + X k )~ li2 is the symmetric positive definite square root of (I-hXj-t 4- 
Xk) 1 Now', applying the Jacobian from (Al) of the Appendix by making the 
transformation 


y } =(i + x 2 + .. + x k y m x x (i + x 2 + + x k y U2 => 

(£112 

dYi = \I + X 2 + ... +X t \ 2 dX x for fixed X 2 , , X k ' (10) 

L\ = \l + X 2 + . +X k )\~ (a2+ +a * +1+Pl+ + fc) 
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* Woftr~i ^r i+ +ai+i+Pi+ 

= \l + X 2 + + x k )\^ ++aM+ ^^ 

.. r ,( a i) r ,( a 2 + +a trt +Pi+-. + Pt) 

^(a, +a 2 + .. + a t+1 +p, + ,. + p t ) ( } 1) 


| 

for gl(ai)> - ,?l(a 2 +. , + a w + .. + P (t ) > The integral is evaluated by using 

a real matrix-variate type-2 beta integral of Appendix (A6) Now, integral over X 2 
yields the following, denoting it by L 2 : 

+ . +x k \ <a2+ +ak+ ' +h+ +m ox 2 

= \I+X 3 + ..+x k \- (a3+ + “* +I+P2+ +te) 

x r j ,(a 2 )r j ,(a3 + . +a t+1 + P 2 + --- + P*) 
r p (a 2 +•• +“lr+l +P 2 + - + Pt) 


I 

for iR(a 2 )>^-^ — ,9t(a 3 +... + a k+l +P 2 +. . + P*) > ^ J = ^ Proceeding like 

this we have the final result as follows • 


c -l _ [A r p( a ;) r p( a 2 -fl + " + P; +-- + Plr) > 

{;=! r p ( a ; +- +a t+ i +P 2 + - + Pt) 


( 13 ) 


From the gamma structure in (13) we can have an interesting result . 
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Theorem 1 For arbitrary t h ...,t k the joint product moment of the determinants of 
Xi, , Xfcis given by the following : 


E\\x x nx 2 r. ix*n 

= (tt EpK 1 + " + tX * +1 th t + / 14) 

*{;= 1 r p (a ; +... + a jt+1 +P ; +.. +p fc -i, - J 


with t 0 =0, for > -5R(a,), j=l,. , k, 5R(/i+. +t ; ) < K(a /+ i+.. +a*+i+p ; +. +p*), j 
=1, k. Note that (2 6) is available from c k x of (14) by changing a 7 to a } + tj,j = 1, . k 
Some interestmg results can be derived from (14). Let t\= 0 = . =tj_\ = t J+ . =t k . 
Then from (14) we get 

E[\ jx,| 0 ] 

_ r p^ a j + r r^ a i + ' + - + a ^i + Py +- + Pt 

^(ctj) rp(a J+ i + . +a jt+1 +p^ + . + P*.) - 

/ T r K +' y ~T i ) r ( a ; + i +■ • +a ^l +P; +• +Pi zh zSl 

m=i r(a ; --^4)r(a ;+1 + .. + a k+l +P ; + .. + P* 

= E(Z *). E{Z\) 

where Z [; are independently distributed type-2 beta random variables with the 
Parameters (ct ; a ;+I + ... + a t+1 +P ; + ... + P* - J fJ-),m = l,2. 

Corollary 1 : When X], . ,.JC k have the joint density in (9) then for a specific j, \Xj\ is 
structurally a product of p independent real type-2 scalar beta random variables and 
the t'j moment of |Aj| is given in (16). 


(15) 

(16) 

(17) 
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Theorem 2 . Let W } - Wj > 0, j = 1, k be independently distributed real matrix 

variate type-2 beta random variables with the parameters (a, a ;+1 + ..+a 4+ i+f3+...+ p k ) s j 
= l, ,,k. Consider the tranformation 


X, =(l+w k ) l/2 .(/+W 2 )' /2 W ] (I+W 2 ) 1/2 (l+w k ) m 
X 2 =(I+W k y n . .(I+JV 3 ) U2 W 2 (I+W 3 ) V2 (I +W k ) U2 


(18) 


X k .xHI+W k ) V2 W k .,{I+W k ) V2 


x k = w k 


where (I+Wj) ]a denotes the symmetric positive definite square root of (J+Wj) for j = 
1, ,.,k Then X 2 ,.. Jf*have the joint density as given m (9). 

Proof Since X 2 is free of W u X 2 is free of W h W 2 , and so on, the Jacobian matrix in 
the transformation (18) is of a tramgular block matrix format. Hence the Jacobian is 
available from the product of the determinants of the diagonal blocks. The j th diagonal 
block behaves like with W J+ h W k j = 1, ,.,k fixed. Then for fixed Wfj,..., 

W k , 


dX, = U + W, 


;+i 


p+i P+1 

2 .\l + W k \—dW r j = 1, 


,k. 


Hence, 


,(*-lXp+l), (i-2Xp+l) 

dX, a ... a dX k = |7 + W k \ a 1 1 + W kA | 2 

(£t!) 

x\l + W 2 \ 2 dW,A...AdW k . 


( 19 ) 
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Since W u . ,W k are assumed to be independently distributed the joint density of 
W h ,W k is the product of matrix-variate type-2 beta densities Denoting the joint 
density of W u ,W k by g(fV i,. ,,W k ) we have 

g(W I} ,W k ) dWi a.. AdW k 

= |n^( a ;> a j + l + - +a <r + l+P; + - +P*)| 

xjlll^ f'^l I + W J f ( “^ aj+1+ + “ i+1+Pj+ +Pi) | dW x A . A dW k 


From the transformation in (18) we have the following • 
I+X k = I + W k 

I+X k +X kA = (I+W k ) m (I+W kA )(I+W k ) m 


i+x k +. +x, = (i+w k y , 2 (i+w kA ) V2 

X. .{I+W 2 ) V2 {I+W x )(l+W 2 ) xn .{I+W k ) x ' 2 


Let us examine, 


W 2 -\^ k \ k 2 \i + x 2 + .+x^| pl \i+x k \ 


iPfc-i 


Ir „ „ l' (a > + +a i-+l + Pl + + Pit) 

x\I+X,+ . +X k \ dX^A.AdX k . 
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= i ^ r 2 i w k i“*~ 2 1 1 + w x r (ai+ +ai+i+p > + 

X I 1 j-(‘ x 2+ +“<:+' + P2+ + PO 


x\\ + W k r ( “* +c * +I+Pt > dW x A. . a dW k 

The nght side multiplied by c k is the product of the densities of W h , W k Hence 
X\, , X k have the density as m (9) 

The converse also holds 

Theorem 3 Let X\, , X k . have the joint density as in (9). Consider the transformation 

in (10) for some real symmetric positive definite p x p matrices W h .. , W k Then 
W\, W k are independently distributed, and further that IT) is a matnx-variate real 
type-2 beta with the parameters (a p a ]+ x +. + a t +i+P/ +. .,+P*), j = 1,. k 

Theorem 3 can also be restated in a different form which looks like fully different 
but it is essentially the same transformation in a different form. 

Theorem 4 ■ Let Xi, , X k be jointly distributed as in (9) Consider the transformation 

W k = X k 

w k . x = (I + x k y U2 x k -d(i + x k y m Y 

w k . 2 = (/ + x*+ x k -.\)~ xn x k -i[(i + x k +x k - x y U2 Y 

: ( 20 ) 

w x = (/ + x k +...+ x 2 y V 2 x x [(i + x t +.. +x 2 y V2 Y 

Then W u .... ; W k are independently distributed where IF) has a real matnx-variate type- 
2 beta density with the parameters (a ; , a, +! +.. + a i+ ]+P ; +.. +Pt), j = \, . ,k. 
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Proof Here we will use the property that a real symmetric positive definite matrix A 
can be written as A = A l ' 2 (A 112 )' where A 1 ' 2 need not be symmetric but square and 
nonsmgular From the transformation in (20) observe the following • 

X k =W k =$I+X k =I+W k 

X k „ x =(I+X k ) 1,2 W k . x [(I+X k ) 1 /2 ] '=(/+ W k ) 1/2 W k _ ] [(I+ W k ) m ] 1 
X k - 2 =(I + x k + w k _ 2 [(i + X k + x k . x ) V2 y 


But 


I+X k +X k _i =1+ w k +(l+ W k ) 1 12 W k _ , [(/+ W k ) 1 /2 ] 1 

-(i+w k ) xn ((i+ w k _\)[(i+w k ) U2 y 
Hi+w k ) 1,2 (i+w k ^) l,2 [(i+w k _ l ) U2 y[(i+w k ) V2 y 


Therefore we may take 


=(I+X k +X k _ , ) 1/2 =(/+ w k ) 1 n (/+ W k _ , ) 1 1 /2 


and 


=l(/+x k +x k . l ) I/2 y = [(/ + w k _j) 1/2 ]' [(/+ w k ) U2 y 


Then 


X k - 2 -(/+ W k ) 1 12 (/+ W k _ , ) 1 n W k _ 2 [(/+ W k - 1 ) 1 ' n ] ' [ (/+ 1 /2 ] ' 

and so on, and finally 

X,={I+W k ) m (I+W k _ x ) X12 ... (/+r 2 ) 1/2 W l [(I+Wi) m ]'...[(I+W k ) m y 
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Thus 

i+x k + +x 2 +x x =(i+w k y /2 ... (i+w 2 ) 1/2 (i+w,)[(i+w 2 ) V2 y [{i+w k ) m y 

But this is the same transformation m (19) and hence by Theorem 3 the result follows 

Note that if W } has a real matnx-vanate type-2 beta distribution with the 
parameters (a ; , a, +i+ . a*+l+p ; + then IF' 1 must have a real matnx-vanate 

type-2 beta distribution with the parameters (a J+) + . +a,t+i+P./ + +p fc a,) But from first 
principles it is quite difficult to prove this result. We will state this as a theorem 

Theorem 5 . Let, (X\, , X k ) have the joint density as in (9) Consider the following 
transformation . 

V k = X~ x 

v k ~i = (/ + x k y hl X~\[{I + x k y m Y 
v k -2 = (/ + x k + x k - x y U2 [(/ + x k +x,_ 1 >" 1/2 ]’ 


( 21 ) 


Vi - (i + x k + + x 2 y V2 x p 1 [(/ + x k + +x 2 y U2 y 


Then V,, ,V k are mutually independently distributed, and further, that V } is a real 
matnx-vanate type-2 beta with the parameters (a^+i +. +ai+i+p ; + . +p to a ; ),. ), j = 1,. ., 
k 

Proof We will provide an independent derivation of this result rather than observing 
that V, is of the form V } = Wf x of Theorem 4 Observe the following : For any real 

symmetric positive definite matrix A one can find a square nonsingular matrix C such 
that A = CC', |C| 0 In the following steps we will denote this C by A m and write A 

= {A V2 )(A m )' From the transformation in (21) we have the following : 
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x k =v;' ^i+x k =i+v-' 

I+V r {I+X k + ..+X ]+I ) V2 (I+X?W+X k +.. +X J+1 ) 1/2 ]' 
for ] = 1, Also 


x k ^=(i + x k y n vMu + x k ) m y 
=(i + v-'y' 2 v- \[(i+v k - l ) U2 y 


x k -2 = (I + x k + x k . x y m v~! 2 [(i + x k +Xi_i) _1/2 ]' 

= (/+ Vj-' ) 1/2 ( 1 + v k ~\ ) 1 /2 [( i + ) 1/2 ] • (i+F t -') ,/2 ]' 

Xi= ((/+ v k ~') U2 ... (i+v; x ) V2 vy x ...[(i+v?) m y ... W+v; x ) xn y 

Further, 

I + X k =I + V k ~ x 


i + x k +x k . x = (/+ f ; 1 ) 1 /2 ( i + v -\ ) [ (/+ f ; 1 ) 1 72 ] 1 


I + x k + ..+X x =(/+ Vy x ) xn ...(\ + V 2 x ) xl 2 (I+V[ x ) 


xtO + F^) 172 ]*... [(l + F t -1 ) m ] 
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Hence, 




p + 1 »+ 1 »+i 

.iF,r ^ iz+^r^ . i/+F 2 r^ 


jH-l 

x*r 2 =m 




. p+i ^ p+i 

' k 


_p+i 


2 ...i v k i 2 i/+F*r 2 i/+F 3 r 2 


P+l n+l 

I ~=| Ft f fc+ ~ 


1 1 +x k i fe_1 =i v. r fe -> i / + k i Pt -' 


| / + X 2 + X fc _, | p *- 2 =| V k _ x |' Pi - 2 1 v. j - p *- 2 \I + V k | p *~ 2 |/ + Ft., | fe - 2 


|/+Z,+. + X 2 | p >=|F k r Pl .|F 2 |- p '|/ + FJ p > ,.|/ + F 2 | Pi 


|/+x t +...+x,| 


“(«l+ +a*+l+Pl+ +p*) 


H Ft | 


“i+ +“i+i+Pi+ +Pt 


oti+ +at+i+P!+ +Pt 


x|i+F f (0C1+ + “* +1+Pl+ +P * J 


xh+f.i 


-(«i+ 


+a t+i+Pi + + po 


The Jacobian matrix is of a triangular format and the Jacobian is the product of the 
determinants of the diagonal block matrices. It can be seen to be the following : 
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dX x a a dX k =\V k 


r(k ^E±L (p+1)rr _ ( *. 2) £±L (i , +1) 


W„ 


k-\ 


.(*-!> 


p + \ 


P +1 


/+F. I 2 .. I/ + FJ 2 dK A A dV k . 


Now, substituting all these we obtain the following 

(p+ i) a ill 

cj^T 2 ..\x k \ k 2 |/+x 2 +..+xj ( V+x 3 + l P2 


I r ,, r v r (ai+ +a ^l +p l + 

X | / + | 1 1 / + X k + . . + | dX^ A. A dX k 


k 

ns 

;= i 


«/+!+ +a*+i+P;+ +pjf 


F, 


ill 

2 


I/+F 


-(a,+ +a* + i+P;+ +0*0 


( 22 ) 


where 

c - n>(a l + " +a t+1 + P; + ••• + Pi-) 

S( a j)n( a ;+i + - +a i+i + Pj + — + P*) 


(23) 


for 9l(a 7 )>ijl,; = l, ,*,9t(a 7+1 +. +a t+I + P, + .p*)>ijl, 7=1, ... k. Tins 

completes the proof that VJs are independently distributed and that each V } is a real 
matrix-vanate type-2 beta matrix random vanable with the parameters 
( a ;+ 1 +"' a i-+l +P; + •■ + Pi><*;)> J- 1, • , k- 

Appendix 

Some of the Jacobians of matrix transformation, integrals over scalar functions of 
matrix argument and type-1, type-2 beta integral will be listed here. For details on 
these as well as for related results see Mathai 1 . 
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Y=AXA', X=X\ \A\ *0=>dY=\A\ p+l dX 


(AI) 


where the nonsmgular matrix A is free of the elements in the real symmetric matrix 
X=X' 

Let T=(t, } ), tj, >0, j= 1 , ,,p, t,j = 0 for ? <y be a lower triangular matrix with 
positive diagonal elements Consider a pxp real symmetric positive definite matrix 
X=X'>0 such that X=TT' Then 


f p 

X = TT' =>dX = 2 p \Y\t JJ 

b=> 


p+i-j 



(A2) 


With the help of (A2) one can evaluate the matnx-vanate gamma integral in (7) and 
obtain the result in (6) Let L be a p x p nonsmgular matrix of functionally 
independent real variables Let X-X' Then 

Y=X~\ X=X', \X\ ^G=^dY=\X\' ip+x) dX (A3) 


With the help of (Al), (A3), (7), (6) we can evaluate the real matrix-variate beta 
integral and obtain the following results ' 


5 P (a,P) 


r p (a)T^(P) 

r p (a + P) 




I-X 


I 2 


dX 


1 o<*=y</ 



p + i 

I-X\~ 2 dX 


/ x=x’>0 



1-X r (a+w dX 


(A4) 


(A5) 
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= 1 


x=x‘>0 



\I-X |- (a+p) dX 


(A6) 


where (A5) and (A6) are called real matrix-variate type-1 and type-2 beta integrals 
respectively 
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Abstract 

The equilibrium of a self gravitating cylmdncal polytrope with a general magnetic field and 

rotation has been discussed 

(Keywords polytropes/ self-gravitating) 

Introduction 

There has been interest in the study of self-gravitating systems in different 
geometrical configurations 1-6 . Recently, study of perfectly conducting, self gravitating 
infinite cylinders (and jets), of compressible as well as incompressible fluids, has 
gained renewed interest, for example by Radwan 6 as it has relevance in describing and 
finding condensations within astronomical bodies especially the spiral arms of the 
galaxies. Most of the workers have assumed a homogeneous density for fluid cylinder. 
In order to generalize this for variable density cylinders we have chosen a polytropic 
equation of state. Even though the effect of geometry can be visualised, explicit 
results are not available except for spherical geometry. A spherical self- gravitating 
polytrope has been discussed in detail by Chandrashekhar 1 , who has observed that the 
equilibrium size of the sphere approaches infinity for n=. 5; n being the polytropic 
index. In the present work, we intend to find such limits when the configuration is a 
cylindrical one and determine the effect of a general magnetic field (with toroidal and 
poloidal components) and rotation on it. Incidentaly it has been found that the 
cylinders do not become boundless even for values far beyond n=5. 

Basic Equations 

Here we consider an ideally conducting polytropic fluid cylinder pervaded by a 
general magnetic field of the form H = (0fJ z (W ) , H$(co) where H$ and H z are the 
toroidal and the axial magnetic field components respectively and co denotes the 
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radial coordinate. The configuration is assumed to be rotating with a uniform angular 


velocity V about the z-axis. 

The equations determining the equilibrium (Sood 3 ) are 


-VP+pVO-oVx(Vx a") + — (V xH) x H = 0 

4 n 

(1) 

V 2 0=-4tcGp 

(2) 

and 


V. H=0 

( 3 ) 


where V=Vt J ,I z is the unit vector along z-axis and ®,p and p are the gravitational 
potential, kinetic pressure and density respectively 

It is convenient to express H as 

H = + (4) 

a cm 


where T and P are toroidal and poloidal scalars 

The Lorentz force (except for a constant multiplier) thus becomes 


_ 1 _ 

4tt 


(VxH)xH = — h 

An 


- a-^(a 2 p)+ 4 ~ (a 2 P) - a T 
aa a d(S 



( 5 ) 


Choosing the equation of state in the form p-K p 1+1/n , p is the pressure and p is the 
density of the fluid and AT is a constant, we obtain [after taking divergence of equation 
( 1 )] 


\_d_ 
co d& 


a^;^(« + i)p ,/ '-} 


+ AnGp - IV 2 - — V.[V x H) x H)] = 0 
4 ti 


( 6 ) 



EQUILIBRIUM OF SELF GRAVITATING POLYTROPES 


355 


If we choose the magnetic field to be purely toroidal and also take 


T d /rT2 


(o' T) = K Vi 


we obtain 


V — (VxH)xH = — AT, 
_4np J 2n 


For a purely axial magnetic field we choose 

— — (//") = — K 2 
2pcte 2 % 

Now assuming that the magnetic field has both the axial and the toroidal 
components, the eqn 6) takes the simple form 

1 A- ^ A. (K{n + l)p 1/B )+ 47t G p - 2F 2 + — (K, + K 2 ) = 0 
cd da da 2 tc 

With p=k@", this equation further reduces to the following form 


-i i f t ' 

-- - E — © +©''’-£=0 

W ^ J 


where ^ is a dimensionless variable defined by / w=a% with 


a = („ + l)X 1/(n+I) and e= -^-~[2F 2 ~(K X +K 2 )1 2%}. The eqn. (8) for the 

_4itG J 47cG 

special case e=0 is 


-—[^—©1 + ©" =0 

J 
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which is similar to the Lane-Emden equation / or spherical polytropes and is to be 
solved under the initial condition © =1 and dldt, © =0 at \ =0 The first condition 
defines X as the density along the axis of the cylinder while the second condition 
follows directly from eqn (9) 


Existence of Solutions 

The eqn (9) is analytically soluble for «=1 yielding the solution ®(^)=7 0 (^) with 
its first zero at = 2 4048. For other values of n, the equation has to be solved 
numerically It is obvious that the value of e plays an important role in the solution of 
the equation 

(i) For 6=0, where the magnetic energy balances the rotational energy, we find 
that the value of n hardly exceeds 30. 

(n) For positive values of e that is when rotational energy dominates the 
solution is found only in a limited region as depicted in Fig 1 

(in) The solution is found to exist for all negative values of e as expected. Here 
magnetic energy is the dominant factor 



Fig 1 - Existence of solution 

However, it is found that the radius of the cylinder hardly undergoes any change with 
increasing polytropic index when e < -0 5 Since the equation becomes non-linear for 
n > 1 , the above results are based on numerical computations even though the results 
for n = 1 can be obtained analytically. 
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Physical Characteristics 
The radius of the cylinder is given by 


ffl. =o£i ={K(n + \) 


1 

4/t G 


] m X°- n) 2n^ 


where is the first zero of ©(^) 
The mass per unit length is 


M- 2np&da =-^ 1 [K(n + l)A. 1/ ” 




and the mean density is given by p = M/(na + 2 ) so that p/X = — (— ©) 5=?1 


Results and Conclusions 

(1) The variation of the radius of the cyhnder( denoted by § 1 ) is shown m Fig. 2 
for e = 0 We can conclude that for n > 25, almost approaches infinity. In 
fact the variation of the size of the cylinder with n can be described by a 
quartic function of the form 

§i = 2.414 - 0.3176 (m-I) + 0 3973 («-l) 2 - 0 04351 («-l) 3 +0.00245 (n-1) 4 
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(u) From the pclytropic equation of state, a large value of n implies a near 
isothermal behaviour. An analytical solution exists only for n = 1 m this 
geometry as well as in spherical geometry For r, - 5 the boundary in the 
spherical case approaches infinity whereas the boundary is finite for the 
cylindrical case and may be regarded infinite for much larger values of r as 
mentioned above This difference may be assigned to the different 
topologies The effects of magnetic field and rotation are depicted explicitly 
by the parameters e < 0 and e > 0 respectively (Fig. 3 ) 



Poiytropic index n 



Radia! distance Fig. 4- Adiabatic case 
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(m) Variation of p with the radial distance i; for different e has been depicted in 
Fig 4 where it is observed that the equilibrium size of the cylinder decreases 
as e decreases i e the magnetic field confines the boundary In Fig 5 effect 
of h on density has been shown graphically for 6=0 It is interesting to note 
that for large values of n nearly the complete mass of the cylinder shrinks 
close to axis even though the boundary of the cylinder is still large (Fig 1) 
bearing a semblence to possible black holes m this geometry Changes m the 
i ahie of e give results on the expected lines as is evident from Fig 5 



0 0.5 t 1.5 2 2.5 3 

Radial distance 

Fig 5- Effect of n on density 
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Abstract 

Transient decrease in cosmic ray intensity followed by a slow recovery is called a Forbush 
decrease event m cosmic ray studies Such two Forbush decrease events were noted first m 
August 19 to 29,1999 another in February 15 to 23, 1999, which occured during high solar 
activity period of solar cycle 23 It is noteworthy that both the events were accompanied with 
several intense solar flares and geomagnetic storms. In this work we have done a systematic 
study to derive relationship between cosmic ray decrease and increase m geomagnetic activity 
We have found a significant relationship between cosmic ray intensity and decreasing profile of 
Dst values Results of analysis indicate a strong relationship between geomagnetic activity and 
cosmic ray intensity on short-term basis 

(Keywords forbush decrease/solar flares/geomagnetic activity) 

Introduction 

A number of cosmic ray researchers established that Forbush decrease events m 
cosmic rays are produced by perturbation in interplanetary conditions and that these 
perturbations originate either from solar flares or from magnetic field structures 
associated with interplanetary solar wind streams 1 ' 2 The perturbations could be 
produced from shock waves, coronal mass ejections and flare generated high speed 
solar wind streams 3, 4 . In 1975 Barouch and Burlaga 5 have reported that the high 
magnetic field regions m interplanetary space are associated with Forbush decreases 
Further it has been demonstrated that these cosmic ray decreases are not related to 
turbulent or random motions in the field, while the regions of high field strength m 
interplanetary space are found to be responsible for causing Forbush decreases. These 
regions consist of magnetic blobs and magnetic clouds ejected from active solar 
regions. Inteiplanetary shocks have comparatively ordered field structure The 
turbulent field in the environment of shocks and high-speed streams are simply 
tangential discontinuties 5 



362 


PANKAJ K SHRIVASTAVA and PRASANTI SHRIVASTAVA 


The idea of spiral cone like region, which extends along the interplanetary 
magnetic field, is also investigated In this study we have studied the characteristics of 
Forbush decrease events of August 19 to 29,1999 and February 15 to 23,1999 m 
relation with geomagnetic field variations Recently, it has been investigated that the 
transient disturbances m interplanetary medium such as magnetic clouds and bi- 
directional events are also producing cosmic ray decreases on short-term basis 6 

Data and Method of Analysis 

Forbush decrease events of August and February 1999 were identified from the 
hourly plots of the cosmic rays intensity We have used the daily values of the two 
neutron monitors Kiel (2 23 GV) and Climax (2.97GV) Solar flares, Dst indexes, Ap 
index, Solar sunspot number were taken from the Prompt reports of Solar Geophysical 
Data books. Data related with geomagnetic activity/solar activity for the event periods 
were also taken from the Solar Geophysical Data book 

Results and Discussion 

Fig. 1 shows the Forbush decrease event of August 1999 Onset time of this event 
is in early hours of August 19, 1999 and it remains low during four to five days. This 
figure shows cosmic ray intensity plots for seven neutron monitors- Calgary, 
GooseBay, Moscow', Climax, Kiel, Beijing and Haleakala covering a wide range of 
cut off rigidities. Several intense solar flares were observed during the event period, 
solar flares of August 17 which occurred in Noaa region of 86658 seem to be more 
powerful m producing decreases in cosmic ray intensity Further solar flare events of 
August 19, 20, and 21 m the Noaa region are also support to enhance the decrease m 
cosmic ray intensity Major solar flare events of August 20 and 21 are responsible to 
produce maximum decrease in August 23. A senes of geomagnetic events of August 
16,17,18 and 22 provide the possible solar terrestrial relationship during the period of 
Forbush decreases 

Similarly Fig 2 shows the plot of Forbush decrease everij: of same neutron 
monitor stations for the period February 15 to 23,1999 Onset time of this event is on 
February 17 and it remains low during two days. A recovery period of this event is 
four days. We noticed two major solar flares on February 15 and 16 on Noaa region of 
8462 which may have affected to produce the large decrease in cosmic ray intensity. 
We have also noted geomagnetic storms on February 17 and 18 which coincide with 
maximum decrease of cosmic ray intensity. 



Percent 



Fig 2 - Plots for hourly values of cosmic rays for several neutron momtors for February 1 999 event 
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Fig 3- Daily values of Sun Spot numbers (Rz) and geomagnetic activity (Ap) for August 1999 



DAYS 

Fig. 4-* Same as Fig 3, but for Feburary, 1999 



Aug. 1999 



Feb.1999 


Fig 5- Variation of Dst values dunng the event periods of August 1999(left) and February 1999(right) 







RELATIONSHIP BETWEEN COSMIC RAYS & GEOMAGNETIC ACTIVITY 


365 


Further we have plotted the daily values of Ap index for both the events as shown 
in Fig (3 & 4) respectively Daily values of Ap for the event of August 1999 are 
plotted on Fig 3 and similarly these geomagnetic indices are plotted in Fig 4 for 
another Forbush decrease of February 1999. High Ap values are seen for both the 
event periods 

Conditions of interplanetary medium play an important role in energy 
transmission from solar wind to magnetosphere that significantly modulate cosmic ray 
particles Forbush decrease events as we observed during August and February 1999 
are examples of such events to draw a possible relationship between cosmic ray 
intensity and interplanetary features. 

We observed high solar wind speed and interplanetary magnetic field B ( nT) to 
event periods Occurrence of several intense solar flares along with flare generated 
solar wind streams are expected as causes of Forbush decrease event of August 1999. 
High values of B ( nT) coincide with the event period and low solar magnetic fields are 
also observed for event period of February 1999 Solar wind speed data are not 
available during event period to draw any meaningful conclusions Hence we can infer 
the same interplanetary conditions for both the events periods. However, variations in 
cosmic ray intensity of event periods depend on the larger/smaller disturbances in 
interplanetary conditions. High-speed solar wind streams particularly solar flares 
generated are found more responsible factor in producing transient decrease in cosmic 
rays This result is consistent with previous results reported in 1996 by Shrivastava 
and Shukla 4 by taking two different types of high-speed solar wind streams. 

To show the relationship between geomagnetic activity and cosmic ray intensity 
we have plotted the daily values of geomagnetic Dst index for both the Forbush 
decrease event periods. Daily values of Dst index are plotted in Fig 5 for the events of 
August 1999 and similarly the values of Dst index for the event period of February 
1 999 are also plotted in right panel of the figure. A significant variation in Dst values 
is seen for both of events periods. Variational profiles of Dst index for these two 
events periods vary according to variational profiles of Forbush decrease events. A 
large decrease m Dst index for the February 15 to 23 coincides with main phase of 
Forbush decrease event of February 1999. Similarly we also observed a significant 
decrease m Dst index which coincides with decrease of cosmic ray intensity for entire 
period of Forbush decrease event of August 1999. In both the events, we observed 
significant decrease in cosmic ray intensity as well as geomagnetic Dst index during 
the main phase of storm. 
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The relationship between the solar flare and cosmic ray intensity decrease 
observed dur in g these two periods are similar to that in previous investigations 7 ’ 8 

One sudden storm commencement (SSC) was observed on February 17,1999 in 
association with the Forbush decrease event of February 1999 As we know SSC is a 
signature of arrival of shocks waves. Therefore it can be inferred that shock waves 
also play an important role in producing large cosmic ray decreases 

Significant decreases were seen in geomagnetic Dst index during both the periods 
of Fobush decreases. The values of Dst index show significant transient decrease m 
similar pattern as that of Forbush decrease. The geomagnetic disturbance index Dst at 
low latitude comes from the outward blowing zonal current system called the ring 
current. It is assumed that the massive compression of the magnetosphere and 
enormous intensification of the large-scale magnetosphric current system reflected in 
Dst leads to a significant effect on cosmic ray measurement near the earth 

Coronal mass ejections (CMEs) are known as one of the causes of interplanetary 
disturbance because they involve a large amount of mass and energy ejected into 
interplanetary medium In a number of recent studies, CMEs are found as a 
responsible factor in cosmic ray modulation and geomagnetic disturbances 9 . We have 
observed a series of CMEs events during the period of both the Forbush decrease 
events We observed two CMEs events on August 17 & 22, 1999 Similarly a series of 
CMEs events were observed on February 12, 13, 14, 20 & 21. Now, it can be inferred 
that the CMEs were also one of the factor in producing decrease in cosmic rays and 
increase in geomagnetic activity. 


Conclusions 

It has been concluded from the analysis that- 

1. Both of forbush decrease events are found in association with increase in 
geomagnetic activity. 

2. The main phase of cosmic ray decrease coincides with maximum Dst values 

3. Solar flares and coronal mass ejections are found to be responsible factors in 
producing decrease in cosmic ray intensity and increase in geomagnetic 
activity. 

4. One can distinctly see these on the basis of sudden storm commencement 
(SSCs). Occurrence of SSCs event of February 17, 1999 clearly exhibits the 
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effect of interplanetary shock waves on cosmic ray intensity which shows a 
step like decrease with fast recovery in shorter time span. 
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